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ABSTRACT 

This  report  presents  the  theoretical  basis  and  validation  measurements  supporting  a 
new  configuration  of  Fourier-transform  IR  spectrometer  optimized  for  differential  detection. 
The  optimized  configuration  exploits  the  optical  subtraction  capabilities  of  a  double-beam 
Michelson  interferometer.  It  has  been  designed  to  satisfy  two  requirements,  the  suppression 
of  the  spectral  background,  which  overlays  the  desired  spectral  signature,  and  the  suppression 
of  the  internal  IR  emission  contributions  arising  from  the  spectrometer  itself.  The  resulting 
configuration  is  a  fully  symmetrical  FTIR  double-beam  interferometer  with  a  special 
beamsplitter  assembly.  To  verify  the  proposed  design  a  prototype  has  been  developed  and 
tested  experimentally.  The  results  obtained  with  this  prototype  indicate  a  high  level  of 
suppression  with  a  residual  usually  smaller  than  2  %  in  the  spectral  range  from  3.5  to  15  pm. 
The  analysis  of  the  results  indicates  that  the  remaining  residual  can  be  explained  by  small 
dissymetries  in  transmittance  and  thickness  of  the  two  ZnSe  substrates  forming  the 
beamsplitter.  Finally,  this  spectrometer  configuration  allows  a  permanent  radiometric 
calibration,  which  represents  a  definitive  advantage  over  conventional  FTIR  instruments. 


RESUME 


Ce  rapport  presente  les  bases  theoriques  et  les  mesures  de  validation  qui  sous-tend  une 
une  nouvelle  configuration  de  spectrom&tre  IR  par  transformation  de  Fourier  (FTIR)  qui  est 
optimisde  pour  la  detection  diffdrentielle.  Cette  configuration  optimisde  exploite  les  capacitds 
de  soustraction  optique  de  l’interfdrometre  de  Michelson  4  double  entrde.  II  a  6t6  congu  de 
manure  b.  satisfaire  deux  impdratifs,  la  suppression  de  i’arriere-plan  spectral  qui  masque  la 
signature  ddsiree  et  la  suppression  de  remission  IR  auto-induite  par  le  spectrom&tre.  La 
configuration  rdsultante  est  un  interfdrom&tre  FTIR  &  double  entree  parfaitement  symdtrique 
muni  d’un  dispositif  particulier  de  lames  sdparatrices.  Afin  d’dvaluer  le  design  propose,  un 
prototype  a  6t6  mis  au  point  et  v6rifi6  expdrimentalement.  Les  rdsultats  obtenus  avec  ce 
prototype  indiquent  un  degrd  dlevd  de  suppression  avec  un  rdsidu  infdrieur  h  2  %  dans  la 
region  spectrale  3,5  &  15  pm.  L ’analyse  des  rdsultats  indique  que  ce  rdsidu  peut  s’expliquer 
par  de  petites  dissymetries  entre  les  transmissions  et  les  dpaisseurs  des  deux  substrats  formant 
le  dispositif  sdparateur.  Finalement,  cette  configuration  de  spectromdtre  permet  un 
etalonnage  radiomdtrique  permanent,  ce  qui  reprdsente  un  avantage  marque  sur  les 
spectrometres  FTIR  conventionnels. 
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EXECUTIVE  SUMMARY 


The  Canadian  Forces  are  using  a  variety  of  electro-optical  systems  for  the  remote 
sensing  of  target  characteristics.  In  the  past  few  years,  the  development  of  passive  IR  spectral 
sensors  has  grown  to  a  point  at  which  it  is  now  considered  a  vital  technology  for  the  remote 
monitoring  of  battlefield  environment,  providing  unique  information  on  the  ongoing 
maneuvers.  For  instance,  one  of  the  most  promising  applications  in  passive  spectral 
surveillance  is  the  remote  detection  and  identification  of  toxic  vapors.  It  has  been  shown  that 
such  a  detection  can  be  performed  remotely  using  a  Fourier  Transform  InfraRed  (FTIR) 
spectrometer  which  probes  the  IR  emission  spectra  of  the  associated  gases.  This  technique  was 
originally  conceived  to  provide  ground  soldiers  with  an  early  warning  of  chemical  threats. 

Defence  Research  Establishment  Valcartier  (DREV)  is  currently  developing  a  passive 
FTIR  technique  for  the  passive  remote  monitoring  of  battlefield  gaseous  emissions.  There  are 
two  known  difficulties  encountered  in  the  application  of  such  a  technique.  First,  it  requires 
the  suppression  of  the  spectral  background  which  overlays  the  desired  spectral  signature. 
Second,  it  requires  the  handling  of  the  internal  IR  emission  arising  from  the  spectrometer 
itself.  To  overcome  these  two  limitations  DREV  has  recently  proposed  a  novel  approach  of 
differential  detection  based  on  a  specially  designed  double-beam  FTIR  spectrometer.  In  this 
system,  two  beams  coming  from  different  scenes  can  be  spatially  combined  on  a  single 
detector  and  optically  subtracted  in  real-time.  Thus  if  one  beam  entering  the  interferometer 
corresponds  to  the  target-plus-background  scene  and  the  other  beam  corresponds  to  the 
surrounding  background,  then  the  resulting  spectrum  is  free  of  background  and  self¬ 
instrument  emission.  This  capability  is  particularly  suited  for  the  detection  of  weak  IR 
emission  signatures  in  a  strong  atmospheric  background. 

This  report  presents  the  theoretical  basis  and  the  validation  measurements  supporting 
the  design  of  the  proposed  double-beam  FTIR  spectrometer  for  differential  detection.  The 
research  effort  has  been  concentrated  in  formulating  the  proper  requirements  for  an 
optimized  configuration  and  finding  practical  solutions.  In  particular,  to  achieve  these 
requirements  it  has  been  proposed  to  use  a  fully  symmetrical  FTIR  double-beam 
interferometer  with  a  special  beamsplitter  assembly.  To  verify  the  proposed  configuration  a 
prototype  has  been  developed  and  tested.  The  results  obtained  with  this  prototype  indicate  a 
high  level  of  suppression  with  an  instrument  residual  usually  smaller  than  2  %  in  the  spectral 
range  from  3.5  to  15  pm.  Moreover,  it  has  been  found  that  this  configuration  permits  a 
permanent  radiometric  calibration,  which  represents  a  definitive  advantage  over  conventional 
FTIR  instrument  for  real-time  sensing. 

The  results  of  this  work  may  also  have  a  significant  impact  on  several  other  remote 
sensing  applications  involving  atmospheric  pollution  monitoring  from  the  ground  or  from  an 
airborne  platform.  In  addition,  this  work  contributes  directly  to  Canadian  inputs  to  TTCP 
JTP-13  and  NATO  RSG-8  international  technical  groups. 
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1.0  INTRODUCTION 

There  exist  a  variety  of  electro-optical  (EO)  sensors  being  developed  to  assist  the 
Canadian  Forces  in  surveillance  operations.  Among  them,  passive  IR  spectral  sensors  appear 
to  become  a  vital  technology  for  the  remote  sensing  of  battlefield  environment  which  can 
provide  unique  information  on  the  ongoing  maneuvers.  For  instance,  spectral  information 
might  be  used  to  remotely  probe  different  types  of  equipment,  camouflages,  smokes  and 
clothing.  One  of  the  most  promising  military  application  in  passive  spectral  surveillance  is 
the  remote  detection  and  identification  of  toxic  vapors.  It  was  originally  conceived  to  provide 
ground  soldiers  with  an  early  warning  of  chemical  threats.  It  has  been  shown  that  such  a 
detection  can  be  performed  remotely  using  a  Fourier  Transform  InfraRed  (FTIR) 
spectrometer  which  probes  the  IR  emission  spectra  of  the  associated  gases  (Refs.  1-3). 
However,  there  are  two  main  difficulties  associated  with  such  a  passive  spectral  technique. 
First,  it  requires  the  suppression  of  the  spectral  background  which  overlays  the  desired 
spectral  signature,  and  second,  it  requires  the  handling  of  the  internal  IR  emission  arising 
from  the  spectrometer  itself.  In  this  report  we  propose  an  approach  to  address  these  two 
limitations. 

Defence  Research  Establishment  Valcartier  (DREV)  is  currently  developing  a  passive 
FTIR  technique  for  the  passive  remote  monitoring  of  battlefield  gaseous  emissions.  The 
DREV  approach  takes  advantage  of  the  differential  detection  provided  by  a  specially 
designed  FTIR  dual-beam  interferometer  with  adjacent  Fields  of  view.  In  this  system  two 
beams  coming  from  different  scenes  can  be  optically  combined  on  a  single  detector  and 
subtracted  in  real-time.  Thus,  if  one  beam  entering  the  interferometer  corresponds  to  the 
target-plus-background  scene  and  the  other  corresponds  to  an  equivalent  surrounding 
background  scene  alone  then  the  resulting  target  spectrum  is  almost  free  of  background. 
This  is  particularly  suited  for  detection  of  weak  IR  emission  signatures  in  a  strong 
atmospheric  background. 

A  major  limitation  encountered  with  standard  FTIR  spectrometers  is  that  the  output 
signal  contains  parasite  radiation,  i.e.  the  internal  IR  emission  arising  from  the  spectrometer 
itself.  Consequently,  a  special  calibration  procedure  involving  two  reference  measurements 
must  be  applied  to  obtain  and  correct  for  this  self-emission  term  in  addition  to  the  usual 
responsivity  parameter.  In  current  FTIR  spectrometers  this  self-emission  term  is  generally 
uncontrolled  and  variable  in  time.  Thus,  the  best  way  to  acquire  radiometrically  accurate 
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spectra  is  by  doing  frequent  calibration  measurements.  This  constraint  of  frequent 
calibration  represents  a  serious  burden  reducing  the  efficiency  and  the  applicability  of  FTIR 
instrument  especially  for  real-time  sensing. 

In  this  report  we  describe  an  FTIR  spectrometer  optimized  for  the  differential 
detection  (Refs.  4-5).  This  special  configuration  permits  the  simultaneous  suppression  of  the 
spectral  background  and  the  instrument  self-emission  spectrum.  In  addition  it  allows  a 
permanent  calibration,  which  represents  a  definitive  advantage  for  real-time  remote  sensing. 
In  Chapter  2,  we  first  review  some  general  concepts  related  to  Fourier  transform  spectrometers 
(FTS).  The  proposed  FTIR  configuration  and  the  related  design  equations  are  presented  in 
Chapter  3.  In  Chapter  4,  the  theory  of  optical  interference  is  used  to  specify  a  beamsplitter 
configuration  that  matches  the  requirements.  In  Chapters  5  and  6,  the  results  obtained  with  a 
prototype  are  presented  and  discussed  as  a  validation  of  the  proposed  configuration.  Finally, 
Chapter  7  summarizes  the  work  and  draws  conclusions  on  the  anticipated  performance  of 
this  optimized  FTIR  spectrometer. 

This  work  was  performed  at  DREV  between  April  1996  and  November  1997  under 
Thrust  2d  -  Land  Forces  Tactical  Surveillance  and  Counter-Surveillance,  Work  Unit  2da25: 
Gaseous  Emission  Monitoring  for  Surveillance:  Feasibility  Study. 

2.0  BACKGROUND  INFORMATION 

2.1  FTS  Generalities 

In  this  section,  we  review  some  basic  aspects  related  to  Fourier  transform  tpectrometers. 
More  information  on  the  subject  can  be  found  for  instance  in  Ref.  6.  The  basic  characteristic 
of  an  FTS  is  to  act  as  a  modulator.  The  principle  is  that  the  interferometer  modulates  the 
incident  beam  intensity  according  to  the  wavelength  of  the  radiation.  The  interferometer 
most  often  used  is  a  Michelson.  A  schematic  of  this  interferometer  is  presented  in  Fig.  1. 
The  beam  from  the  source  is  divided  by  a  beamsplitter  (BS)  into  two  separate  components 
which  are  recombined  after  traveling  different  optical  path  lengths.  In  a  Fourier 
spectrometer,  the  Michelson  interferometer  is  designed  in  a  way  that  one  or  both  of  its 
mirrors  can  move  periodically  (note  that  in  Fig.  1,  only  one  mirror  is  moving)  around  a  mean 
position  creating  a  periodic  optical  path  difference  between  the  two  interfering  beams.  In  this 
case  when  a  monochromatic  radiation  of  unit  amplitude  (Fig.  1)  is  incident  on  the 
interferometer  the  emergent  amplitude  A(x)  is  given  by 


UNCLASSIFIED 

3 


A(x)  =  (r  t)  +  (t  r)  e  1  2nvx  [  1  ] 

where  v  is  the  wavenumber  (inverse  of  the  wavelength,  1/X),  x  is  the  optical  path  difference 
between  the  two  interfering  beams  and  r  and  t  represent  the  amplitude  reflection  and 
transmission  coefficients  of  the  beamsplitter,  respectively.  Note  that  the  temporal  variation  of 
the  amplitude  terms  in  eq.  1  has  been  neglected  assuming  an  observation  time  long  enough  to 
take  the  averaged  amplitude.  This  simplified  radiation  source  (monochromatic  radiation  of 
unit  amplitude)  is  helpful  in  establishing  the  basic  characteristic  of  the  instrument,  i.e.  the 
optical  responsivity  K.  For  the  moment  we  neglect  the  electronic  responsivity.  As  seen  below, 
the  responsivity  K  depends  essentially  on  the  transmission  and  reflection  of  the  many  optical 
elements  encountered  in  the  optical  path.  By  definition,  the  emerging  intensity  I(x)  is  found 
by  multiplication  of  the  amplitude  A(x)  with  its  complex  conjugate  A(x)*  yielding 

I(x)  =  A(x)  A*(x)  =  2RT  +  RT  (e  *  2nvx  +  e  ->  2 rcvx)  r2] 


FIGURE  1  -  Principle  of  the  Michelson  interferometer  used  as  a  modulator  in  a  Fourier 
transform  spectrometer 
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where  the  beamsplitter  reflection  and  transmission  in  intensity,  R  and  T  are  given  by  (r  r*) 
and  (t  t*),  respectively.  The  interferogram  F(x)  is  defined  as  the  modulated  component  of 
the  above  intensity,  i.e. 

F(x)  =  K  (e  1  2nvx  +  e  -i  2nvx)  [3] 

and  where  K  =  RT  has  been  defined  as  the  optical  responsivity.  Equation  3  represents  the 
interferogram  obtained  from  a  monochromatic  source  passing  through  an  ideal 
interferometer. 


When  a  polychromatic  beam  of  intensity  B(v)  is  incident  on  a  FT  interferometer  it 
generates  an  interferogram  that  can  be  represented  by  the  continuous  summation  of 
individual  modulations  associated  to  each  spectral  element,  i.e. 


F(x) 


B(v)  K(v)  (e 


i  2ttvx 


+  e  -i  2twx)  dv  , 


[4] 


where  the  responsivity  K(v)  is  now  wavenumber  dependent.  Using  a  mathematical  artifice  (it 
does  not  affect  the  physical  representation)  which  consists  in  extending  the  domain  of 
integration  to  include  the  negative  values  of  v  and  assuming  symmetric  spectra  for  B(v)  and 
K(v)  where 


B(-v)  =  B(v)t. 

[5] 

1 

<_ 

III 

2 

[6] 

we  find,  after  rearrangement,  an  interferogram  expression  more  appropriate  for  processing. 
In  this  case  and  by  application  of  the  basic  Fourier  transform  property  we  have 


F(x) 


■i: 


B(v)  K(v)  e  «  2nvx  dv 


and 


[7] 
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B(v)  K(v)  =  j  F(x)  e  -i  2*vx  dx  . 

J-co  [g] 

Equation  7  states  that  the  interferogram  F(x)  generated  by  an  ideal  instrument  is  given  by  the 
inverse  Fourier  transform  of  the  raw  spectrum  B(v)  K(v)  while  eq.  8  states  that  the  raw 
spectrum  is  given  by  the  Fourier  transform  of  the  interferogram. 


As  opposed  to  an  ideal  instrument  a  real  interferometer  usually  introduces  an  additional 
phase  shift  \|/(v)  due  to  an  imperfect  compensation  from  dispersive  elements  (e.g.  substrates 
of  beamsplitter  and  compensator).  In  this  case  the  pair  of  Fourier  transform  becomes 


F(x)  =J  [B(v)  K(v)  e  i  V(v)]  e  * 2nvx  dv 


and 


[B(v)  K(v)  e  i  V(v)]  =  j  F(x)  e  -i  2twx  dx  . 


[9] 


[10] 


Equation  10  is  the  basic  relation  used  in  Fourier  spectroscopy  where  the  complex  raw 
spectrum  is  simply  given  by  the  complex  Fourier  transform  of  a  real  double-sided 
interferogram.  Note  that  the  absolute  radiometric  spectrum  B(v)  is  not  a  direct  output  and 
the  instrument  responsivity  must  be  calibrated  from  auxiliary  measurements  performed  on 
standard  reference  sources.  In  today’s  instrument  the  computation  of  complex  spectra  based 
on  eq.  10  is  not  a  burden  anymore.  Fast  Fourier  Transform  algorithm  (FFT)  coupled  to 
actual  capabilities  of  small  PC  computer  can  deliver  a  high  resolution  spectrum  almost  in  real 
time. 


In  summary,  it  is  important  to  emphasize  that  a  real  instrument  is  characterized  by  its 
complex  responsivity 

K(v)  =  K(v)  e  ■  V(v)  [n] 

where  K(v)  and  \j/(v)  represent  the  module  and  the  phase  of  the  instrument  responsivity  K(v). 
The  bold  face  notation  is  used  throughout  this  paper  to  identify  complex  numbers.  The 
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remaining  part  of  this  document  is  mainly  dedicated  to  the  analysis  of  the  instrument 
responsivity  for  different  instrument  configurations.  In  particular,  the  inter-relation  between 
self  IR  emission,  instrument  responsivity  and  calibration  is  treated  in  view  to  find  optimum 
configurations.  In  this  context,  a  useful  relation  to  compare  one  interferometer  configuration 
with  another  is  found  by  inserting  eq.  1 1  into  eq.  4,  which  gives 

f(x)  =  K(v)  cos  [  2ttvx  +  y(v)  ]  [12] 

where  f(x)  is  the  interferogram  component  corresponding  to  a  single  spectral  element  of  a 
spectrally  uniform  source  of  unit  intensity  defined  by  a  8-Dirac  i.e.,  B(v)  =  5  (v*— v)  This  is 
the  relation  used  here  as  a  reference  for  comparing  the  impact  of  internal  component 
configurations  on  the  interferometer  responsivity. 

2.2  Self-emission  and  Calibration  of  Standard  FUR  Spectrometers 

The  schematic  of  the  instrument  shown  in  Fig.  1  corresponds  to  a  standard  single-input 
beam  Michelson  interferometer.  For  a  real  instrument  of  this  type  the  output  signal  is 
actually  composed  of  radiation  from  different  origins.  As  it  turns  out,  an  important  part  of 
this  radiation  comes  from  the  external  source  under  study  L.  The  remaining  contributions 
are  generated  by  the  instrument  itself.  This  parasite  radiation  is  due  to  graybody  emissions 
and  stray  light  reflections  from  the  interferometer  components  (lenses,  mirrors,  inner  walls 
and  others).  It  is  referred  to  as  the  self-emission  of  the  instrument.  Figure  2  identifies  the 
various  sources  of  radiation  intervening  in  a  single  measurement:  L  is  the  spectral  radiance  of 
the  external  source  (target  source),  SEjn  and  SEout  represent  the  self-emissions  (raw  spectra) 
generated  by  the  input  and  the  output  optics,  respectively,  and  SEgS  represents  the 
beamsplitter  self-emission  (raw  spectrum).  In  the  infrared  the  self-emission  terms  mainly 
depend  on  the  temperature  of  the  intervening  components.  It  is  quite  significant  especially 
for  instrument  operated  at  or  near  room  temperature.  Taken  separately  each  of  these 
radiation  sources  generates  a  specific  interferogram  and  consequently  a  corresponding  raw 
spectrum.  The  sum  of  these  individual  raw  spectra  is  given  by 

S  =  (K  L)  +  SEin  +  SE0Ut  +SEBs  [13] 

which  can  be  simplified  by  defining  an  instrument  offset  O  such  that 
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FIGURE  2  -Source  and  self-emission  contributions  from  interferometer  components.  L  is 
the  spectral  radiance  of  the  external  source,  SEjn  and  SE0ut  represent  the  self¬ 
emission  contributions  from  the  input  optics  and  the  output  optics,  respectively, 
and  SEbs  represents  the  beamsplitter  self-emission  contribution. 


0=X(SEin  +  SEout  +  SEBs), 

k  [14] 


yielding 


S  =  K(L  +  0)  [15] 

Equation  15  represents  the  usual  way  to  link  the  complex  raw  spectrum  S  to  the  target 
spectrum  L  (the  one  to  be  evaluated)  taking  into  account  the  instrument  characteristics  which 
are  the  complex  spectral  responsivity  K  and  the  complex  spectral  offset  O.  As  seen  in  eq.  14 
the  spectral  offset  of  the  instrument  is  a  rather  complicated  expression  involving  three 
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different  self-emission  terms  that  are  theoretically  difficult  to  evaluate  particularly  for  this 
type  of  instrument  (single-beam). 

Usually  the  self-emission  terms  are  not  stable  in  time  due  to  thermal  variations  of 
instrument  components.  For  this  reason,  the  best  way  to  acquire  radiometrically  accurate 
target  spectra  is  by  doing  frequent  calibration  measurements.  Ideally,  this  calibration  is 
performed  concurrently  with  each  source  measurement.  The  current  method  used  for  the 
accurate  calibration  (radiometric)  of  interferometric  spectra  is  known  as  the  two-temperature 
calibration  method.  The  best  variant  of  this  method  was  proposed  by  Revercomb  et  al.  (Ref. 
7).  In  this  method,  two  reference  blackbodies  of  known  radiances  Bhot  an(I  Bamb  are 

required  to  solve  the  two  unknowns,  the  spectral  responsivity  and  the  spectral  offset  of  the 
instrument.  Often  in  FTIR  work,  the  hot  reference  blackbody  is  kept  at  a  constant 
temperature  near  60°  C  and  the  second  one  is  at  ambient  temperature  near,  20°  C.  Defining 
the  two  corresponding  raw  spectra  of  reference  as  Shot  and  Samb  then  the  solution  of  eq.  15 

yields 


K  = 


Shot  ~  Samb 

Bhot  “  Bamb 


[16] 


and 


O 


ibBhot  -  ShotBamb 

Shot  "  Sarnb 


[17] 


And  the  calibrated  target  spectrum  is  obtain  from 


L  = 


[18] 


The  above  discussion  illustrates  the  usual  limitation  (radiometric  calibration)  encountered 
with  most  of  standard  FTIR  instruments,  i.e.  the  presence  of  variable  and  uncontrolled  self¬ 
emissions,  which  necessitates  three  consecutive  measurements  (target  and  two  references)  for 
obtaining  a  single  calibrated  spectrum  L.  This  represents  a  burden  especially  for  real-time 
applications.  In  order  to  make  the  instrument  more  stable  radiometrically  and,  consequently, 
more  simple  to  calibrate  we  propose  to  exploit  and  optimize  certain  attributes  of  the  double¬ 
input  beam  interferometer  design.  As  seen  below  this  type  of  interferometer  is  also  suited  for 
developing  an  FTIR  spectrometer  adapted  for  differential  detection. 
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FIGURE  3  -Schematic  view  of  a  double-beam  Michelson  interferometer  (a)  that  uses  comer 
reflectors  rather  than  standard  flat  plates.  Figures  b  and  c  illustrate  the  two  out-of- 
phase  interferograms  intervening  in  a  double-beam  interferometer  (see  text). 
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2.3  Double-input  Beam  Interferometer 

Among  the  few  numbers  of  double-beam  interferometers,  the  modified  Michelson  with 
comer  cube  retroreflectors  represents  a  good  design  to  start  with  from  the  point  of  view  of  its 
optimization:  self-emission  suppression  and  simplified  calibration.  This  is  the  basic  design 
retained  in  our  approach.  Figure  3a  represents  a  schematic  view  of  the  double-beam 
Michelson  interferometer.  The  use  of  comer  reflectors  rather  than  standard  flat  plate  mirrors 
has  two  main  advantages  for  our  application.  First,  the  reflected  beams  from  comer  reflectors 
are  offsetted  laterally  from  the  incident  beams.  This  has  the  great  advantage  of  decoupling 
input  and  output  channels,  as  opposed  to  the  standard  single-beam  interferometer  (Fig.  2) 
where  incident  and  reflected  beams  are  superimposed.  This  attribute  makes  it  possible  to 
adjust  independently  the  self-emission  of  each  input  ports  without  affecting  in  any  way  the 
output  optics,  which  is  not  possible  with  the  single-beam  configuration.  Note  that  the 
instrument  also  has  a  second  output  optics.  For  the  sake  of  simplicity,  the  following  analysis 
is  restricted  to  only  one  of  the  two  outputs  but  our  conclusions  apply  equally  to  both  in  a 
straightforward  manner.  With  this  type  of  interferometer  the  two  beams  from  the  two  inputs 
are  optically  combined  in  a  subtractive  manner.  This  optical  subtraction  arises  from  the  fact 
that  the  two  interferograms  associated  with  each  input  port  can  be  made  out  of  phase, 
depending  on  the  instrument  design  as  discussed  in  the  following  section.  Figures  3b  and  3c 
illustrate  the  two  out  of  phase  interferograms  intervening  in  this  double-beam  interferometer. 

3.0  THE  PROPOSED  FTIR  SPECTROMETER  CONFIGURATION 

The  strategy  behind  the  proposed  FTIR  spectrometer  consists  in  designing  a  double-beam 
interferometer  such  that  the  two  self-emission  terms  from  each  input  ports  fully  cancels  out 
by  optical  subtraction,  i.e.  when  the  two  self-emission  terms  are  equal  in  intensity  and  180° 
out  of  phase  leading  to  a  theoretically  perfect  cancellation.  To  our  knowledge,  such  a  design 
to  suppress  the  self-emission  in  FTIR  spectrometer  by  optical  subtraction  has  not  been 
proposed  in  the  past.  This  constitutes  the  basic  effort  of  the  present  work:  an  FTIR 
spectrometer  design  which  is  optimized  to  suppress  the  internal  self-emission  by  optical 
subtraction  using  a  properly  balanced  dual-beam  interferometer. 

The  proposed  design  is  the  result  of  a  previous  investigation  on  the  subtraction 
performance  of  a  commercial  dual-beam  interferometer,  the  MB- 100  manufactured  by 
Bomem.  Although  the  MB-100  was  not  specifically  optimized  for  optical  subtraction,  its 
utilization  provided  a  good  understanding  of  the  basic  principles.  From  that  investigation  it 
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appeared  that  there  are  two  main  factors  limiting  the  accuracy  of  the  subtraction.  The  first 
limiting  factor  is  related  to  the  optical  configuration  of  the  two  input  ports  (mirrors  and 
windows).  These  two  inputs  must  be  as  symmetric  as  possible  to  ensure  the  same  beam 
attenuation  and  consequently  the  same  emissions  from  both  ports.  For  instance  if  the 
absorption  of  a  mirror  (or  lens)  from  input- 1  is  higher  than  its  counterpart  from  input-2,  then 
the  resulting  unbalance  corresponds  to  a  larger  emission  from  input- 1.  In  practice,  it  is 
relatively  easy  to  satisfy  this  symmetry  requirement  with  a  proper  mechanical  design  and  the 
use  of  identical  optical  components  in  each  input  port. 

The  second  limiting  factor  is  related  to  the  IR  absorption  and,  consequently,  the  self¬ 
emission  of  the  beamsplitter.  Of  the  two  limiting  factors  the  beamsplitter  one  is  the  most 
important.  For  this  reason,  the  remaining  part  of  this  report  mainly  focus  on  the  beamsplitter 
issue  assuming  that  the  two  input  ports  of  the  double-beam  interferometer  are  symmetric. 

More  formally,  the  following  equations  serve  to  establish  the  expected  attributes  of  the 
optimized  instrument.  For  an  FTIR  spectrometer  based  on  a  double-input  port  interferometer 
the  total  signal  S  represents  the  sum  of  two  individual  complex  raw  spectra,  E]  and  E2, 
associated  with  each  inputs.  If  input- 1  and  input-2  aim  at  target  sources  of  radiance  Li  and 
L2,  respectively,  the  resulting  raw  spectrum  is  given  by  (see  Fig.  3a) 

S  =  Ei  +  E2  +  Egs 

with 

Ei=  KiOU  +  SEjm) 

E2  =  K2(L2  +  SEin2) 

where  K],  K2  and  SEini,  SEjn2  are  defined  as  the  responsivities  and  the  self-emissions 
associated  with  input- 1  and  input-2,  respectively,  and  Ebs  represents  the  self-emission  from 
the  beam  splitter.  Note  that  by  comparison  to  the  source  radiance  terms  L]  and  L2,  the  self¬ 
emission  terms  may  have  a  phase  delay  because  they  may  originate  from  different  location. 
To  study  the  radiometric  symmetry  or  balance  of  the  dual-beam  interferometer  it  is 
advantageous  to  proceed  by  analogy  with  the  single-beam  calibration  approach  (eqs.  16-18). 
In  the  particular  case  where  the  calibration  is  performed  on  input- 1  while  source  L2  is  fixed 

for  the  three  measurements  (target  and  two  references)  then  the  corresponding  raw  spectrum 
can  be  written  as 


[19] 


[20] 

[21] 
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S  =  K!(Li  +  Oi) 


[22] 


with 

0,=  (sEinl+|lL2  +  ^SEin2+^).  [M] 

In  this  case  the  source  L2  acts  a  part  of  the  instrument  offset.  Equations  22  and  23  link  the 
complex  raw  spectrum  S  to  the  target  spectrum  Li  (the  one  to  be  evaluated)  taking  into 

account  the  characteristics  of  the  double-beam  interferometer  which  are  the  complex  spectral 
responsivity  Kq  and  the  complex  spectral  offset  Oi  associated  to  channel- 1. 

Alternatively,  if  the  calibration  is  performed  on  input-2  while  the  source  Li  is  fixed  for 
the  three  measurements,  then  the  corresponding  raw  spectrum  can  be  written  as 


s  =  k2(l2  +  o2) 


[24] 


with 


In  this  case  source  L]  acts  as  a  part  of  the  instrument  offset.  Equations  24  and  25  link  the 
complex  raw  spectrum  S  to  the  target  spectrum  L2,  taking  into  account  the  characteristics  of 
the  double-beam  interferometer  which  are  the  complex  spectral  responsivity  K2  and  the 
complex  spectral  offset  O2  associated  to  channel -2. 

From  equations  23  and  25  it  can  be  seen  that  the  spectral  offsets  Oi  and  O2  which 
represent  the  self-emission  contributions  associated  with  channel- 1  and  channel-2, 
respectively,  can  be  minimized  by  the  fulfillment  of  the  three  following  conditions: 

CONDITION- 1:  Balanced  responsivities, 

K  j  =  -  K2  =  K2  e  i  K  [26] 

CONDITION-2:  Balanced  emissions, 
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SEjni  —  SEjn2 


[27] 


CONDITION-3:  Beamsplitter  transparency 
Ebs  =  0 


[28] 


By  application  of  these  three  conditions,  the  resulting  raw  spectrum  derived  from  eqs.  22  to 
25  reduces  to 

S  =  K  i  (Li  -  L2)  .  [29] 


Equations  26  to  28  define  in  a  mathematical  sense  the  conditions  for  an  optimized  FTIR 
spectrometer  as  proposed  in  the  present  work.  As  seen,  eq.  29  does  not  contain  any  self¬ 
emission  terms  from  the  instrument  itself.  The  resulting  raw  spectrum  S  is  directly 
proportional  to  the  difference  between  the  two  source  radiances  Li  -  L2.  Consequently,  for 

single-beam  operations  this  instrument  requires  the  use  of  a  fixed  reference  source  of  known 
radiance  (Lref  =  L2)  attached  to  the  second  input  port.  This  reference  can  be  a  simple  flat 
plate  blackbody  (Fig.  4)  having  a  thermocouple  for  monitoring  the  plate  temperature  (Tref), 
which  allows  the  computation  of  the  Planck  radiance  L(Tref).  In  this  case,  the  calibrated 
spectrum  Lj  is  directly  given  by 


Li  = 


[30] 


There  are  several  advantages  associated  with  this  approach.  First,  it  considerably 
simplifies  the  radiometric  calibration  where  in  principle  only  one  reference  measurement  is 
needed  to  evaluate  solely  the  spectral  responsivity  Kq  of  the  instrument  (no  offset). 

Moreover,  because  the  responsivity  can  be  correlated  empirically  to  the  beamsplitter 
temperature  (Tbs)  then  the  calibrated  spectrum  can  be  generated  in  real  time  directly  from 

the  monitored  beamsplitter  temperature  and  using  pre-stored  responsivity  coefficients  to 
generate  Kq(TBS)  in  real  time.  In  these  conditions,  the  calibrated  radiance  becomes 


Li  = 


Ki(Tbs) 


+  L(Tref)  . 

1 


[31] 
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FIGURE  4  -  Optical  configuration  for  single  beam  measurements  with  the  optimized  FTIR 
spectrometer:  a  flat  plate 

The  net  result  of  the  proposed  design  is  an  FTIR  spectrometer  that  requires  only  the 
monitoring  of  instrument  temperatures  (beamsplitter  and  flat  plate)  to  generate,  in  real  time,  a 
calibrated  spectrum  according  to  eq.  31.  This  represents  a  definite  advantage  over  the  existing 
FTIR  spectrometers  that  use  three  spectral  measurements  (source  and  two  references)  to 
generate  the  calibrated  spectrum.  In  practical  design,  the  casing  of  the  instrument  can  be 
made  to  integrate  the  flat  plate  and  the  beamsplitter  onto  the  same  mount  such  that  Tbs  = 
Tref  resulting  in  only  one  temperature  monitoring  to  obtain  the  calibrated  spectrum,  i.e. 

Ll  =  (  §  ■■■- +  L(Tbs)  1 . 

\  Ki(Tbs)  )  [32] 

A  final  advantage  of  the  proposed  design  is  for  differential  detection  (using  both  inputs 
as  shown  in  Fig.  3a  )  of  weak  IR  signatures  against  large  spectral  background.  Because  the 
subtraction  is  instantaneously  performed  optically,  the  resulting  background-free  spectrum 
benefits  from  the  full  range  of  signal  amplification  as  opposed  to  sequential  subtraction 
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technique  where  the  range  of  signal  amplification  is  mainly  limited  by  the  background  itself. 
In  this  case,  the  calibrated  radiance  differential  can  be  evaluated  in  real  time  using 


5L  = 


_§ _ 

KjCTbs) 


[33] 


where  5L  =  Li  -  L2-  More  information  concerning  the  differential  detection  of  chemical 
vapors  using  a  double-beam  FTIR  interferometer  can  be  found  in  a  recent  paper  (Ref.  5). 


The  conditions  discussed  above  to  obtain  an  optimized  FTIR  spectrometer  impose  two 
basic  requirements  on  the  instrument  design.  The  first  requirement  implies  that  the 
transmissions,  the  absorptions  and  the  corresponding  self-emissions  produced  by  the 
components  of  the  two  input  ports,  including  mirrors,  lens,  attenuators  and  masks  must  be 
symmetrically  the  same  to  satisfy  CONDITIONS- 1  and  2  (eq.  26-27).  This  can  be  achieved 
in  a  straightforward  manner  using  simple  symmetry  considerations  such  that  one  half  of  the 
instrument  (and  input  components)  can  be  considered  as  the  mirror  reflection  of  the  other 
half  taking  the  beamsplitter  plan  as  reference.  It  is  important  to  note  that  the  thermal 
properties  of  the  instrument  must  also  satisfy  the  same  symmetry  arguments  to  ensure  that  the 
temperature  and  consequently  the  graybody  emission  of  intervening  components  of  each 
input  port  are  symmetrically  the  same. 


The  second  requirement  is  related  to  the  beamsplitter  itself.  In  this  case,  the  design  of  a 
beamsplitter  that  satisfies  the  three  conditions  of  optimization  (see  eqs.  26-28)  cannot  be 
stated  simply.  In  the  next  chapter  we  present  an  analytical  development  which  provides  a  tool 
for  designing  the  beamsplitter.  We  also  propose  an  optimized  beamsplitter  configuration. 

4.0  BEAMSPLITTER  CONFIGURATIONS:  ANALYTICAL  DEVELOPMENT 


In  the  following,  the  basic  theory  of  optical  interference  is  used  to  specify  the  attributes  of 
different  beamsplitter  configurations.  This  analytical  development  begins  with  a  review  on 
some  current  beamsplitter  configurations  that  do  not  match  the  optimization  conditions  (eqs. 
26-28)  followed  by  the  treatment  of  two  cases  having  a  suitable  beamsplitter  configurations. 
The  first  case  corresponds  to  a  beamsplitter  made  of  an  optically  thin  layer  of  a  material. 
Self-supported  thin  films  of  gold,  germanium  and  diamond  are  among  this  category.  It  is 
shown  that  this  configuration  (thin  layer)  adequately  balances  the  two  responsivities,  but  it  has 
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the  drawback  of  generating  unwanted  self-emissions  when  the  layer  is  absorbing.  The  second 
case  corresponds  to  a  beamsplitter  made  of  an  optically  thin  layer  of  air  squeezed  between 
two  thick  substrates  properly  covered  with  antireflection  coatings  on  their  external  faces.  It 
has  been  shown  that  this  beamsplitter  configuration  can  fully  satisfy  the  three  conditions  for 
achieving  an  optimized  FTIR  spectrometer:  balanced  responsivities,  balanced  emissions  and 
layer  transparency.  This  is  the  configuration  proposed  in  the  present  work  as  a  claim  for  the 
optimized  FTIR  spectrometer.  In  the  following,  the  basic  theory  of  optical  interference  is 
used  to  specify  the  attributes  of  different  beamsplitter  configurations.  This  analytical 
development  begins  with  a  review  on  some  current  beamsplitter  configurations  that  do  not 
match  the  optimization  conditions  (eqs.  26-28). 

4.1  Unsuitable  Beamsplitter  Configurations 

One  of  the  most  simple  beamsplitter  assembly  used  in  interferometers  is  made  of  a  single 
substrate  with  or  without  coatings  on  its  interfaces.  Figure  5  illustrates  the  ray  tracing  inside  a 
double-beam  interferometer  that  makes  use  of  such  a  beamsplitter.  It  is  easily  seen  that  the 
interfering  rays  from  input- 1  (Fig.  5a)  propagate  only  once  into  the  substrate,  while  the 
interfering  rays  from  input-2  (Fig.  5b)  propagate  twice.  Consequently,  this  beamsplitter 
configuration  generates  a  permanent  unbalance  equal  to  Ts  (substrate  transmittance)  between 

the  two  responsivities  associated  to  each  input.  In  this  case,  the  only  way  to  satisfy  the 
condition- 1  of  optimization  (eq.  26)  is  by  restricting  the  spectral  band  of  operation  to  a 
region  where  the  substrate  material  is  perfectly  transparent.  This  represents  a  major  constraint 
because  most  of  the  usual  substrates  used  in  infrared,  such  as  germanium  and  zinc  selenide, 
exhibit  absorption  in  certain  parts  of  the  spectrum.  This  limitation  makes  this  type  of 
beamsplitter  unsuitable  for  the  design  of  an  optimized  FTIR  spectrometer.  Another  single 
substrate  configuration  developed  specifically  for  double-beam  interferometers  with  corner 
cube  reflectors  is  shown  in  Fig.  6.  The  main  difference  with  the  above  mentioned 
configuration  is  that  the  number  of  interaction  (transmission  and  reflection)  of  interfering 
rays  with  the  substrate  coatings  is  reduced,  which  minimizes  losses  and  increases  the  sensitivity 
of  the  instrument.  However  the  ray  tracings  show  the  same  limitation  as  for  the  preceding 
configuration.  Interfering  rays  from  input-1  (Fig.  6a)  propagate  only  once  into  the  substrate 
while  interfering  rays  from  input-2  (Fig.  6b)  propagate  twice  generating  the  same  permanent 
unbalance  (Ts)  which  is  not  suitable  for  the  design  of  an  optimized  FTIR  spectrometer. 
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FIGURE  5  -  Ray  tracing  for  a  beamsplitter  made  of  a  coated  substrate:  a-  the  interfering 
rays  from  input- 1  propagate  only  once  into  the  substrate,  and  b-  the 
interfering  rays  from  input-2  propagate  twice  into  the  substrate 
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2 

FIGURE  7  -  Ray  tracing  for  a  beamsplitter  assembly  made  of  a  single-side  coated  substrate 
with  a  compensating  plate  (substrate),  a-  A  part  of  the  radiation  from  input- 1 
exhibits  a  substrate-layer  reflection  (r1)  and,  b-  for  input-2  the  radiation 
exhibits  only  air-layer  reflections  (r). 


UNCLASSIFIED 

20 

A  more  symmetrical  beamsplitter  configuration  is  shown  in  Fig.  7  where  a  splitting 
coating  is  present  on  the  inner  face  of  a  first  substrate  and  a  second  one  is  added  as  a 
compensating  plate.  In  this  case,  the  number  of  passes  through  the  substrates  is  identical  for 
both  rays  from  input- 1  and  input-2.  However,  this  configuration  exhibits  another  source  of 
dissymmetry  which  comes  from  the  coating  itself.  As  seen  by  comparing  Figs.  7a  and  7b,  a 
part  of  the  radiation  from  input- 1  exhibits  a  substrate-layer  reflection  (r'),  while  for  input-2 
the  radiation  exhibits  only  air-layer  reflections  (r).  It  can  be  shown,  from  optical  thin  film 
theory,  that  the  two  reflection  coefficients  r  and  r'  are  not  identical  when  the  coating  is 
absorbing.  This  has  for  consequence  to  unbalance  the  responsivities  and  the  self-emissions 
between  the  two  inputs  in  spectral  regions  where  the  splitting  coating  is  absorbing. 
Unfortunately,  most  of  the  coatings  actually  used  for  IR  beamsplitters  may  absorb  in  certain 
infrared  regions.  For  instance,  germanium  and  thorium  fluoride  coatings  absorb  in  the  far  IR. 
Consequently,  this  beamsplitter  configuration  and,  by  extension,  any  other  configurations 
which  do  not  match  the  number  of  r  and  r'  reflections  between  input-1  and  input-2  do  not 
appear  suited  for  an  optimized  FTIR  spectrometer.  This  suggests  the  use  of  a  fully 
symmetrical  beamsplitter  in  which  r'  =  r  and  where  the  number  of  passes  through  the 
substrates  are  equivalently  the  same  for  input- 1  and  input-2.  This  is  the  strategy  proposed  in 
the  next  sections. 

4.2  Thin  Layer  Beamsplitter 

The  theory  of  optical  interferences  is  used  to  establish  the  responsivities  Ki  and  K2  of  a 

double-beam  interferometer-spectrometer  having  a  thin  layer  beamsplitter.  In  this  case,  the 
ray  tracings  in  Figs.  8a  and  8b  indicate  that  the  output  amplitudes  (Aj  and  A2)  for  beams  of 

unit  amplitudes  incident  on  input- 1  and  input-2  are  given  by 


A 1  =  r  t  +  t  r  e^  , 

[34] 

A2  =  r  r  +  1 1  e^  , 

[35] 

where  r  and  t  represent  the  amplitude  reflection  and  transmission  of  the  layer  and 


<j)  =  2jtvx 


[36] 
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expresses  the  phase  difference  between  the  two  interfering  beams  in  terms  of  the  optical  path 
difference  x  (see  Fig.  8a).  Following  Heavens  (Ref.  8),  the  amplitude  reflection  and 
transmission  of  the  layer  can  be  explicitly  evaluated  by  summing  the  contributions  from 
multiple  reflections  and  transmissions  inside  the  layer.  For  a  layer  symmetrically  bounded  by 
identical  media,  such  as  in  Fig.  8c,  the  summations  yield 


r  _  r1  ( 1  -  e-2i5) 

-  [37] 

1  -  if  e-2i5 

[38] 

1  -  if  e-2i5 

with 

8  =  27rvnd  cos  (0) 

[39] 

where  8  is  the  phase  shift  due  to  a  layer  of  optical  thickness  nd,  n  is  the  refractive  index  of  the 
layer  and  0  is  the  angle  of  refraction  inside  the  material..  Equations  37  and  38  are  obtained 

using  the  following  identities  (Heavens,  Ref. 

associated  to  the  interfaces  of  the  layer 

8,  p.  51-  57)  between  Fresnel  coefficients 

r'l  =  -  rj 

[40] 

t’t  ti  =  (l  -  rf) . 

[41] 

Equating  the  denominators  of  eqs.  37  and  38 

leads  to  a  simple  expression  connecting  the 

amplitude  reflection  and  transmission  such  that 

•- 

; 71 

t  =  r  H  e  2 

[42] 

where  H  is  an  intermediate  variable  defined  as 

t- » 

1 

r— <i 

li 

X 

-  [43] 

2  rj  sin(S  ) 
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FIGURE  8  -Ray  tracings  showing  the  output  amplitudes  (Aj  and  A2)  for  beams  of  unit 
amplitude  incident  on  input- 1  (a)  and  input-2  (b)  and  the  amplitude 
reflection  and  transmission  (c)  of  a  thin  layer  beamsplitter 
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Equations  42  and  43  are  useful  relations  which  simplify  the  analysis  for  obtaining  the 
responsivities  of  the  interferometer. 

Coming  back  to  eqs.  34  and  35  and  multiplying  the  two  output  amplitudes  by  their 
complex  conjugates  (A*i  and  A *2)  yields  the  two  output  intensities  Ii  and  I2  associated 

with  each  input  port,  i.e. 


1 1  =  A 1  A*i  =  (rt  +  trei(!>)  (r*  t*  +  t*  r*  e  -  i  $) , 

[44] 

I2  =  A 2  A*2  =  (  r2  +  t2  e  i  0)  (  r*2  +  t*2  e  - '  <t>) . 

[45] 

Taking  the  coefficient  of  reflection  in  intensity  of  the  layer  R  =  r  r*  and  inserting  eq.  42 
into  the  two  preceding  equations  yield  for  the  two  output  intensities 

II  =(2R2  H2  +  2R2H2  cos  (<j>) ) 

[46] 

I2  =  (  R2  +  R2  H2  +  2  R2  H2  cos  (<j>  -  jc)  ) . 

[47] 

In  eqs.  46  and  47  the  modulated  components  of  intensities  represent 
associated  with  each  input  port.  It  follows  that  the  two  responsivities  K 1 

with  each  inputs  are  given  by 

the  interferograms 
and  K2  associated 

K  ]  =  2  R2  H2  ,  [48] 

K2  =  -  2  R2  H2  .  [49] 

which  satisfies  the  CONDITION  -  1  (eq.  26)  of  balanced  responsivities  where 

K]=-K2.  [50] 

Although  this  simple  configuration  satisfies  the  balanced  responsivities  condition,  it  is  shown 
below  that  such  a  self-supported  layer  can  generate  unwanted  self-emission  when  the  layer  is 
absorbing  (Section  4.5). 
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4.3  Proposed  BS  Configuration:  Responsivity  Evaluations 

This  configuration  corresponds  to  a  beamsplitter  (BS)  made  of  an  optically  thin  layer  of 
air  squeezed  between  two  thick  substrates  properly  covered  with  antireflection  coatings  on 
their  external  faces.  Again,  the  goal  is  to  establish  the  responsivities  Ki  and  K.2  of  the 

double-beam  interferometer-spectrometer  with  this  beamsplitter  configuration.  For  the 
moment,  we  will  assume  that  the  two  substrates  and  their  antireflection  coatings  may  be 
optically  different  to  simulate  a  possible  mismatch.  The  ray  tracings  in  Figs.  9a  and  9b 
indicate  that  the  output  amplitudes  (A\  and  A2)  for  beams  of  unit  amplitudes  incident  on 

input- 1  and  input-2  are  given  by 


A 1  =  ( r  t  tC2  tS2 )  +  ( tci  tsi  t  r  t^2 1|2  e 1  ^  )  [51] 

and 

A  2  =  (  tc2  ts2  f2  tel  tsl  )  +  (  tc2  ts2  r2  e  1  $  )  [52] 

where  r  and  t  represent  the  amplitude  reflection  and  transmission  of  the  thin  air  layer  and  4* 
(eq.  36)  is  the  phase  difference  between  the  two  interfering  beams.  Variables  tel,  tc2  2nd  tsl, 
ts2  are  the  amplitude  transmissions  (Fig.  9a)  of  the  two  antireflection  coatings  and  the  two 
substrates,  respectively.  To  simplify  the  analysis  it  is  convenient  to  express  these  amplitudes  in 
terms  of  their  respective  modules  and  phases  yielding 


tsl  =  tsi  e  1 
ts2  =  ts2  e  > 
tcl  =  tci  e  >ri 
tC2  =  tc2  e  1 r* 


[53] 


Using  these  expressions  into  eqs.  51  and  52  leads  to  after  rearrangements 

Ai  =  tc]  ts]  r  t  ts2  tc2  [t h  til  ei(2oi  +  2r,  )  +  t^ tl2 e  i ( 2°2  +  2r2 lei*!5]  [54] 


and 
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A2  =  t22  [t2  t|i  t&  e  i  (  2oi  +  2n  )  +  r2  t|2  t^2  e  i  (  202  +  2r2  )  e  i  0  ]  .  [55] 

The  two  output  intensities  Ii  and  I2  associated  with  each  input  port  are  found  by  multiplying 
the  two  output  amplitudes  by  their  complex  conjugates  (A*i  and  A  *2),  which  gives 


11  =  Ai  A*1  =  TC]TS1  R2  Ts2Tc2  [  H2  T&Ti h  +  H2  T hXh  +  H2  TslTclTs2Tc2  cos  (  $  +  xj/  )] 

[56] 

and 

12  =  A 2  A*2  =  R2  Th  Th  [  H4  T h  ?h  +ThTh  +  H2  Tsl  Tcl  Ts2  Tc2  cos  (  $  +  V  -  tc)  ] 

[57] 

where  a  phase  offset  due  to  an  imperfect  compensation  of  the  interferometer  has  been 
defined  as 


y  =  2(a2  +  r2-a1  -  H) 


[58] 


and  other  quantities  in  capital  letters  are  defined  as  follows 


R  =  r  r* 

Tsl  = tsl  t*si 
Tcl  =  tcl  t*cl 
TS2  =  ts2  t*S2 
TC2=  tC2  t*c2 


Reflection  in  intensity  of  the  thin  layer  of  air. 

Transmission  in  intensity  of  substrate- 1, 

Transmission  in  intensity  of  the  antireflection  coating  onto  substrate-1, 
Transmission  in  intensity  of  substrate-2, 

Transmission  in  intensity  of  the  antireflection  coating  onto  substrate-2, 


and  H  is  defined  by  eq.  43. 
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FIGURE  9  -  Ray  tracings  showing  the  output  amplitudes  (Ai  and  A2)  for  beams  of  unit 
amplitudes  incident  on  input- 1  (a)  and  input-2  (b)  for  a  beamsplitter  made  of 
an  optically  thin  layer  of  air  squeezed  between  two  substrates  properly 
covered  with  antireflection  coatings  on  their  external  faces 
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Note  that  in  eqs.  56  and  57  the  transmission  in  intensity  of  the  thin  layer  of  air  does  not 
appear  explicitly  but  is  defined  as  T=  R  H^,  which  is  consistent  with  eq.  42  connecting  the 
reflection  and  transmission  coefficients  of  a  thin  layer. 

Finally,  in  eqs.  56  and  57  the  modulated  components  of  the  intensities  represent  the 
interferograms  associated  with  each  input  ports.  It  follows  that  the  two  complex  responsivities 
K  i  and  K.2  associated  with  each  input  are  given  by 


K  !  =K  !  eiV  =  +R2H2T21  T&  T^T^eiV 


[59] 


and 


K2  =  K2  e  iv  =  -  R2  h2  Ts1  Tc1  Tj2  T32  e 

[60] 

which  satisfies  the  balanced  responsivities  CONDITION  -  1  (eq.  26) 
substrates  and  their  antireflection  coatings  are  optically  identical  i.e. 

only  when  the  two 

if 

Fsi  “  Ts2  , 

[61] 

Tci  ~  Tc2  > 

[62] 

then 


K  j  =  -  K  2  =  R2  H2  Tfi  T&  e  i  V  .  [63] 

As  will  be  seen  further,  the  equations  developed  in  this  section  will  appear  useful  in  the 
analysis  of  experimental  results  obtained  with  a  prototype  instrument  designed  with  this  type 
of  beamsplitter. 

4.4  Proposed  BS  configuration:  Emission  of  substrates 

The  goal  of  this  section  is  to  show  that  the  two  emission  components  from  each  substrate 
composing  the  proposed  beamsplitter  configuration  cancels  out  .  First,  it  is  assumed  that  both 
substrates  may  be  absorbing  and  that  their  temperatures  are  identical  yielding  a 
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corresponding  blackbody  radiance  equal  to  Bs.  Using  the  parameters  defined  in  Figs.  9a  and 
10,  it  is  defined  that 


£l  =  (l  -TCJT,i)  Bs  [64] 

where  ej  represents  the  emission  component  from  substrate- 1  which  is  incident  on  the  air 

layer.  Other  quantities  has  been  defined  in  the  previous  section.  Note  that,  in  eq.  64,  if  both 
transmissions  in  intensity  Tci  (AR  coating)  and  Tsi  (bulk  substrate)  are  equal  to  one  then  the 

emission  term  vanishes.  Following  a  ray  tracing  and  a  mathematical  development  similar  to 
the  one  found  in  the  previous  section,  it  can  be  shown  that  the  modulated  component  of  the 
interferogram  associated  with  this  substrate  emission  is  given  by 

MODULATION  (e  i)=  [  (1  -  Tci  Tsl)  Bs  ]  R2  H2  Tsl  Tcl  T22  cos  (4>  +  y) 

[65a] 

or  equivalently  in  the  spectral  domain,  the  corresponding  raw  spectrum  associated  with 
substrate- 1  emission  is  given  by 

RAW  SPEC (e  i)  =  Ki (1  _TclT.sI1  Bs  [65b] 

tcl  tsl 

where  the  complex  responsivity  K]  defined  at  eq.  59  has  been  introduced  in  the  previous 
expression. 

Similarly  the  emission  component  from  substrate-2,  which  is  incident  on  the  air  layer  is 
defined  as 

e  2  =  0  -  TC2  Ts2)  Bs  [66] 

and  the  modulated  component  corresponding  to  substrate-2  emission  is 

MODULATION  (e  2)  =  [  U  -  Tc2  Ts2)  Bs  ]  R2  H2  Tsl  Tci  Th  T*2  cos  (<j>  +  \|/  -  rc) 

[67a] 


or  equivalently  the  corresponding  raw  spectrum  emission  is  given  by 
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FIGURE  10-  Ray  tracings  representing  the  two  internal  emission  components  £i  and  £2 
(see  text)  from  substrate- 1  (a)  substrate-2  (b),  respectively 
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RAW  SPEC  (£  2)  =  K2  (1  ^Tc^T---  Bs  [67b] 

Tc2  1  s2 

where  the  complex  responsivity  K2  defined  by  eq.  60  has  been  introduced.  Inspection  of 
eqs.  65a  and  67a  shows  that  the  two  modulation  terms  are  180  degree  out  of  phase  and 
cancels  out  when  the  two  substrates  and  their  corresponding  AR  coatings  are  optically 
identical,  i.e.  if 

Tci  =  Tc2  t68] 

Tsi  =  Ts2  [69] 

then 

MODULATION  (e  j )  =  -  MODULATION  (e  2) ,  [70] 

which  means  that  the  real  substrate  emissions  are  not  detected  at  the  interferogram  level  when 
the  two  substrates  (and  AR  coatings)  are  optically  identical. 

4.5  Self-Emission  from  the  Splitting  Layer 

The  last  important  source  of  self-emission  arises  from  the  splitting  layer  itself  where 
graybody  radiation  and  stray  light  can  emanate  from  the  layer  and  interfere  with  itself.  In  the 
proposed  beamsplitter  configuration  (Section  4.3)  the  splitting  layer  is  constituted  by  a  thin 
film  of  air  which  is  totally  transparent.  Consequently,  there  is  no  self-emission  expected  from 
this  air  layer.  However,  for  the  sake  of  completeness  we  present  a  simple  analytical  method 
for  evaluating  this  self-emission  contribution  when  the  splitting  layer  is  absorbing.  To  our 
knowledge,  such  a  method  has  not  been  reported  in  the  past.  Our  approach  can  be  used  to 
compute  the  self-emissions  of  a  wide  range  of  splitting  layer  configurations. 

We  will  assume  that  the  splitting  layer  is  absorbing  with  a  complex  refractive  index  n  1  = 
ni  -  iki,  a  thickness  d,  and  is  symmetrically  bounded  by  identical  transparent  media  of  index 
no.  Figures  11a  and  lib  represent  the  layer  together  with  parameters  used  to  model  the 
phenomenon.  The  bulk  material  composing  the  layer  generates  (emission  and  stray  light) 
isotropic  IR  radiation.  Only  the  radiation  inside  the  field  of  view  of  the  detector  will 
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contribute  to  the  interferogram.  The  only  two  components  that  can  contribute  to  the 
interferogram  are  indicated  in  Figs.  11a  and  lib  where  ein  is  defined  as  the  emitted  intensity 

(radiance  units)  of  the  bulk  material.  For  each  direction,  this  emitted  intensity  generates  two 
beams  of  amplitudes  ei  and  e2  as  a  result  of  the  multiple  reflections  inside  the  layer.  After 

reflections  on  the  two  comer  reflectors  such  as  for  usual  interferometer  ray  tracings  (see  Figs. 
8a  and  8b),  the  two  beams  are  recombined  for  interference  at  the  detector  level. 

For  the  emitted  intensity  contribution  in  direction  of  the  fixed  comer  reflector  (Fig.  11a) 
the  amplitude  ei  of  the  first  beam  can  be  evaluated  by  adding  the  multiple  reflections 

intervening  inside  the  layer.  Following  Heavens  (Ref.  8)  and  according  to  the  Fresnel 
coefficients,  defined  in  Fig.  1 1,  this  amplitude  is  given  by 

ei  =  [t'i  +  r'f  t'i  e-2a  e-2iA  +  r'j  t\  e-4a  e-4iA  +...)  [71] 


where 


A  =  2  K  v  n|ff  d  cos  (0i)  [72] 

and 

a  =  2  7t  v  kjff  d  cos  (0i).  [73] 

The  effective  refractive  index  (neff,)  and  effective  absorption  coefficient  (keffj)  can  be 
derived  from  the  Snell  law.  Knowing  that  r’l  =  -  ri,  and  adding  terms  in  eq.  71  yields 


e  l  =  Vei^  7 — TiJ-l0  e'iA 


(l-r? 


»-2a  p-2iA| 


[74] 


which  can  be  expressed  in  a  more  compact  form  using  the  amplitude  transmittance  (t) 
expression  of  the  layer  defined  by  eqs.  38  and  41,  i.e. 


e  i  =  Vein  e“  eiA  . 
ft 


[75] 


From  similar  considerations,  the  amplitude  e2  of  the  second  beam  (Fig.  11a)  is  found  to  be 


^2  =  Ve^  (  r'f  e_a  e~iA  t'i  +  r'f  t'i  e_3a  e_3iA  +...) 


[76] 
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FIGURE  11-  Ray  tracings  for  the  modelization  of  the  self-emission  from  the  splitting  layer. 

Graybody  radiation  and  stray  light  can  emanate  from  the  layer  and  interfere 
with  itself  to  create  a  significant  contribution  to  the  raw  spectrum  (see  text). 
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resulting  in  after  addition 


^2  = 


-  ri  t'i  e~a  e~iA 
(l  -  x\  e-2a  e‘2iA) 


[77] 


which  can  be  expressed  in  a  compact  form  using  the  layer  transmittance  (t)  defined  by  eqs. 
38  and  41, 


[78] 


Consequently,  the  output  amplitudes  corresponding  to  the  two  interfering  beams  associated 
with  this  self-emission  component  (Fig.  1 1  a)  are 


A  i  =  e  i  t  =  Ve hT  —  ea  eiA 


[79] 


A2  =  e2  rei<!)  =  Ye^  tre^ 


which  yields 


[80] 


A  i  +  A  2  =  -M  t  ea  eiA  -  q  re  $) 
ti 


[81] 


The  associated  intensity  la  is  found  by  multiplying  eq.  81  by  its  complex  conjugate  leading 
to 


Ia=  £in^[  T  e2a  +  RRi  -n  ea(  r  t*  ei(<f-A)  +  r*  t  e-i^^))]  [82] 

Ti 

where  capital  letters  representing  intensities  are  defined  as  T  =  1 t*,  R  =  r  r*  and  Ri  =  ri  rj. 

For  the  case  in  which  the  emitted  intensity  contribution  is  in  the  direction  of  the  moving 
comer  reflector  (Fig.  lib),  the  two  amplitudes  ei  and  e2  are  simply  the  reverse  situation  of 

Fig.  11a.  In  this  case,  we  have 


[83] 
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e2  =  Vein  —  ea  eiA  [84] 

ti 

which  yields  for  the  two  output  amplitudes  interfering  at  the  detector  level 

A,=e,t  =  V5-(±L)t2  [85] 

A2  =  e2  r  =  Vej^ ^  ea  e>0  .  [86] 

tj 

Adding  these  two  output  amplitudes  and  multiplying  the  result  by  its  complex  conjugate 
gives  the  corresponding  output  intensity  lb  (Fig.  lib)  incident  on  the  detector 

Ai  +  A2  =  V£h7— ( -n  t  +  r  eaei(‘t> +A) )  [87] 

ti 


and 


Ib  =  £in-:L[Re2a  +  TR1-rie«(  rt*ei('5)+A)  +  r*te-i^  +  A))]  .  [88] 

Ti 


The  quantities  of  interest  in  eqs.  82  and  88  are  the  two  modulated  components  which  add 
together  to  produce  the  interferogram  corresponding  to  the  layer  self-emission.  Defining  this 
self-emission  as  Eiay  it  follows  that 

Ejay  =  ein (-ri  ea) [  r  t* e ^ (e_iA  +  elA)  +  r*te -4> (e_iA  +  e [89] 
Ti 

or 

Elav  =  Bin  ea  2  cos  A  )  [  r  t*  e  ^  +  r*  t  e  40  ]  .  [90] 

Ti 

To  further  reduce  eq.  90  it  is  convenient  to  use  again  the  connecting  relation  between  t  and  r 
(eq.  42)  but  for  and  absorbing  layer  i.e. 

t  =  r  H  ein /2  PI] 


where  now  H  is  a  complex  number  given  by 
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H  =  H  ei  6h  =  Il£l  - _1 - .  [92] 

i  (ea  e_iA  -  e_a  e'^) 

Note  that  eq.  92  is  directly  obtained  by  inserting  the  complex  phase  associated  with  an 
absorbing  layer,  8  =  A  -  ia,  into  eq.  43  and  where  A  an  a  are  defined  by  eqs.  72  and  73, 
respectively.  Taking  eq.  92  into  account  and  using  eq.  91  into  eq.  90  yields 

Eiay  =  ein^r-(-ri  e«2  cos  A  )  [  r  r*  H  e1^  ~J'Qh)  +  r*  r  H  e  '1(<t>  'J '0H)] 

[93] 


or  in  a  simpler  form 

Elay  =  ein  (  ri  ea  2  cos  A  )  (2  R  H)  cos  (  §  +%-  -  Oh  )  [94] 

where  again  R  =  r  r*  and  H  together  with  0H  are  defined  by  eq.  92. 

As  defined,  eq.  94  gives  an  analytical  expression  for  the  self-emission  of  the  layer 
composing  the  beamsplitter.  It  should  be  emphasized  that  eq.  94  is  not  rigorously  exact 
mainly  due  to  the  fact  the  phase  shift  associated  to  the  air-layer  reflection,  ri,  was  neglected  in 

the  development.  This  phase  shift  becomes  important  for  strongly  absorbing  layer.  Thus  eq. 
94  constitutes  a  good  approximation  for  slightly  absorbing  layer  materiel  which  is  usually  the 
case  because  this  is  the  condition  that  maximizes  the  sensitivity  of  the  interferometer.  In 
addition,  for  slightly  absorbing  layer  material  0h  becomes  negligible.  Consequently,  the  self¬ 
emission  interferogram  (  eq.  94)  appears  to  have  an  extra  phase  of  90  degrees  over  the  source 
interferogram. 

Finally,  the  driving  term  in  eq.  94  is  the  emitted  intensity  ejn  defined  as  the  graybody 
radiation  and  stray  light  emanating  from  the  bulk  of  the  material  composing  the  layer.  For  a 
well-designed  beamsplitter  the  scattering  inside  the  splitting  layer  is  usually  negligible 
eliminating  the  stray  light  contribution.  Consequently,  the  only  significant  contribution  to 
the  emitted  intensity  £jn  arises  from  the  graybody  emission  of  the  layer  material,  which  is 

defined  as 
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ein=  (l  -  e'2a)  B(  TBs  )  [95] 

where  B(Tbs)  is  the  Planck  radiance  corresponding  to  the  temperature  of  the  beamsplitter 
Tbs  a°d  2a  is  the  optical  depth  (see  eq.  73)  associated  with  the  layer  material.  In  the 
proposed  beamsplitter  configuration  (Section  4.3),  the  splitting  layer  is  constituted  by  a  film 
of  air  which  is  totally  transparent.  Consequently,  for  the  proposed  beamsplitter  there  is  no 
self-emission  associated  to  the  layer  itself  because 

a  =  0,  [96] 

£in  =  0,  [97] 

and  consequently, 

Elay  =  0  .  [98] 

In  conclusion,  for  the  proposed  beamsplitter  design  the  self-emission  from  the  splitting  layer 
vanishes. 


5.0  EXPERIMENTAL  PROTOTYPE:  RESULT  AND  DISCUSSIONS 


In  order  to  evaluate  the  performance  of  an  optimized  FTIR  spectrometer  i.e.  for  self¬ 
emission  suppression  and  simplified  radiometric  calibration,  we  have  specified  an  instrument 
referred  to  as  the  Compact  ATmospheric  Sounding  Interferometer,  which  was  built  under  a 
contract  with  BOMEM.  The  primary  but  not  exclusive  application  of  this  instrument  is  for 
the  passive  remote  sensing  of  target  vapors  (see  Ref.  5).  In  the  following  section,  a 
description  of  the  CATSI  prototype  and  a  discussion  on  its  measured  characteristics  are 
presented  as  a  validation  of  the  proposed  FTIR  configuration  described  in  this  document. 

5.1  Description  of  the  CATSI  prototype 

The  specifications  of  the  instrument  were  made  to  agree  with  the  three  conditions  of 
optimization  defined  by  eqs.  26  to  28.  For  that,  the  optical  configuration  of  the  two  input 
ports  (mirrors  and  windows)  was  made  as  symmetric  as  possible  to  ensure  the  same  beam 
attenuations  and  the  same  self-emissions  in  both  ports.  The  selected  beamsplitter 
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configuration  was  the  one  presented  in  Section  4.3  made  of  an  optically  thin  layer  of  air 
squeezed  between  two  thick  substrates  properly  covered  with  anti  reflection  coatings  on  their 
external  faces. 

More  precisely,  the  CATSI  prototype  is  made  of  two  identical  4-in.  diameter  Newtonian 
telescopes  optically  coupled  to  the  dual-beam  interferometer.  Figure  12  summarizes  the 
optical  design  and  shows  the  instrument  mounted  on  a  tripod.  An  important  effort  was  done 
to  respect  symmetry  criteria.  The  CATSI  system  allows  measurements  of  spectra  according  to 
the  following  specifications:  scene  fields  of  view  from  4  to  1 1  mrad,  spectral  coverage  from  3 
to  18  jam,  and  a  spectral  resolution  of  1  cm-1  or  greater.  A  flat  plate  mirror  placed  in  front 
of  each  telescope  can  be  rotated  to  the  selected  scene.  The  pointing  capability  of  this  scene 
mirror  allows  azimuth  measurements  from  0  to  180  degrees.  Coarse  adjustments  in  azimuth 
and  elevation  are  simply  achieved  by  rotating  the  whole  assembly  mounted  on  a  tripod.  After 
reflection  on  the  scene  mirror  the  input  beam  is  focused  by  the  Newtonian  telescope  at  the 
entrance  of  the  interferometer  and  then  reflected  by  an  off-axis  parabolic  mirror  to  produce  a 
collimated  beam  of  proper  diameter  in  front  of  the  beamsplitter.  A  double  pendulum 
scanning  mechanism  controls  the  periodic  displacement  of  the  two  comer  reflectors  (CC)  that 
generates  the  interferogram.  The  beamsplitter  consists  of  a  thin  air  gap  (A/4  at  7  pm) 
squeezed  between  two  ZnSe  substrates  having  antireflection  coatings  on  their  external  faces. 
Of  the  two  output  channels  only  one  is  used  at  this  moment.  This  output  module  contains 
parabolic  and  condensing  mirrors  that  focus  the  beam  onto  a  sandwich  MCT-InSb  detector  (1 
mm)  mounted  on  a  microcooler  (EG&G  Judson):  The  MCT  element  is  optimized  for  the  6- 
18  pm  spectral  region,  while  the  InSb  element  is  optimized  for  the  2-5  pm  region.  Two  CCD 
cameras  mounted  on  the  top  of  the  two  telescope  modules  can  be  used  to  aim  and  view  the 
scenes  under  consideration.  With  this  instrument,  two  scenes  coming  from  adjacent  FOVs  are 
optically  combined  at  the  detector  level  yielding  the  spectral  residual  of  the  scenes.  The 
overall  system  forms  a  cube  of  approximately  30  cm  3  and  weights  18  kg. 

5.2.  Results:  Instrument  responsivities  and  residual 

In  Chapter  3  we  have  defined  two  basic  requirements  to  obtain  an  optimized  FTIR 
spectrometer.  The  first  requirement  implies  that  the  self-emissions  produced  by  the  optical 
components  of  the  two  input  ports  (excluding  beamsplitter  components)  must  be 
symmetrically  the  same.  After  a  careful  experimental  investigation  with  the  CATSI  prototype 
we  came  to  the  conclusion  that  there  is  no  major  difficulty  to  achieve  the  desired  symmetry 
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FIGURE  12-  Optical  diagram  and  a  picture  of  the  Compact  Atmospheric  Sounding 
Interferometer  (CATSI) 


between  the  two  input  optics,  i.e.  the  two  related  self-emissions  are  quite  balanced.  Rather,  the 
main  limitation  appeared  to  come  essentially  from  the  beamsplitter  itself.  The  remaining  part 
of  this  document  emphasizes  this  point  by  showing  the  impact  of  the  beamsplitter  accuracy 
on  the  performance  of  the  CATSI  prototype. 

In  Chapter  3.0  it  has  been  established  that  the  raw  spectrum  S  delivered  by  a  double-beam 
interferometer  of  this  type  can  be  defined  referring  to  input- 1  characteristics, 

S  =  K,(L1+Oi),  [99] 

Oi  =  (SEinl+|iL2t|2-SEin2+^-)  [100) 


or  alternately  with  the  characteristics  of  input-2 
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S  =  K2  (L2  +  02),  [101] 

0HSEi"2+£tLi+^SEtai+Td  [102) 

where  we  remember  that  L]  and  L2  are  the  two  target  source  radiances  with  Kq,  K2  and 
SEjnl,  SEjn2  defined  as  the  responsivities  and  the  self-emissions  associated  with  input- 1  and 
input-2,  respectively.  For  this  particular  configuration  of  beamsplitter  that  uses  a  thin  layer  of 
air  as  a  splitting  layer  it  has  been  shown  (Section  4.5)  that  in  this  case  the  self-emission  term 
Elay  =  0  yielding  Ebs  =  0  and  consequently, 


o,  =(sEini+|iL2  +  |lSEin2), 

[103] 

02  =  (SEin2  +  |lL1+KlSEinl}. 

[104] 

Note  that  the  two  substrate  self-emission  contributions  are  included  in  the  input  port  emission, 
SEinl  and  SEin2- 


The  experimental  measurement  of  Kq,  K2  and  O],  O2  provides  a  quantitative  approach 
to  evaluate  the  compensation  characteristics  of  the  interferometer.  For  that  the  two- 
temperature  calibration  method  described  in  Chapter  3.0  was  applied.  Figures  13a  represents 
the  modules  of  the  measured  responsivities  for  input- 1  and  input-2.  The  good  agreement 
between  the  two  spectra  from  700  to  3000  cm"^  confirms  the  high  level  of  symmetry 
obtained  with  the  proposed  FTIR  spectrometer  configuration.  Fig.  13b  compared  the 
modules  of  two  measured  offsets  obtained  when  the  calibration  is  done  on  input-1  (Oi)  or 
alternatively  when  it  is  done  on  input-2  (O2)  for  the  case  in  which  both  inputs  look  exactly  at 
the  same  source  of  radiance  L  =  Li  =  L2  (blackbody  at  15°  C).  From  inspection  of  eqs.  103 

and  104  it  is  concluded  that  this  good  agreement  between  the  two  offsets  and  L  in  the  spectral 
region  from  700  to  3000  cm"'  is  consistent  with  these  conditions  Kq  =  -  K2  and  SEinl  = 
SEin2  yielding  O]  =  -  L  and  O2  =  -  L,  as  observed  in  Fig.  13b. 

Another  way  to  look  at  the  optical  subtraction  capability  of  the  instrument  is  provided  by 
a  suitable  recast  of  eq.  101  (taking  input-2  as  a  reference)  in  a  more  convenient  format,  i.e. 
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S=  K2(5L  +  re2), 


[105] 


where  5L  =  L2  -  Li,  and  re2  accounts  for  a  possible  optical  asymmetry  and  is  defined  as  the 
instrument  residual  given  by 


rH,+tfe)MSEi”2+^SEi4 


[106] 


In  this  case,  the  calibration  equation  for  a  dual  beam  interferometer  operating  in  differential 
detection  mode  is  given  by 


6L'(i^-r4 


[107] 


Figure  13c  shows  an  example  of  the  measured  instrument  residual  re2  (bottom  curve). 
For  an  ideal  instrument  this  residual  should  be  zero  throughout  the  spectrum.  For  our 
prototype  it  has  been  found  that  the  instrument  residual  (unbalance)  is  not  exactly  zero  but 
usually  smaller  than  2  %  in  the  region  from  700  to  3000  cm"^  (3.5  -  15  p.m).  This  is 
approximately  a  factor  of  10  smaller  than  the  instrument  residual  measured  on  a  conventional 
double-beam  interferometer  (BOMEM  MB-100)  in  certain  spectral  region.  However,  there  is 
an  unexpected  strong  residual  in  the  from  500  to  700  cm-*  (see  Fig.  13c).  After  an  extensive 
investigation  we  came  to  the  conclusion  that  these  imperfections  (2  %  and  strong  residual) 
come  from  the  limited  performance  of  the  manufactured  beamsplitter  (prototype).  An 
undesired  dissymmetry  between  the  optical  properties  of  the  two  ZnSe  substrates  forming  the 
beamsplitter  was  found.  This  is  partly  due  to  non-identical  antireflection  coatings  present  on 
the  external  faces  of  each  substrate,  and  partly  due  to  a  difference  between  the  transmission 
associated  with  each  substrate.  Both  effects  introduce  a  dissymmetry  in  the  self-emission  and 
the  attenuation  associated  with  each  input.  A  useful  model  which  predicts  the  instrument 
residual  re2  in  terms  of  the  substrate  characteristics  is  presented  and  validated  in  Section  5.3. 

In  this  model  the  instrument  residual  is  directly  given  by 


rej“(1+i^)[Ll'B‘]' 


[108] 
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This  relation  entirely  defines  the  instrument  residual  in  terms  of  the  beam  splitter 
temperature  (Bs),  the  source  radiance  (Li)  and  the  responsivity  ratio.  Although  the  ultimate 

performance  of  the  CATSI  prototype  has  not  been  reached  yet  due  to  manufacturing 
imperfections  in  the  actual  beamsplitter,  the  previous  relations  indicate  that  these 
imperfections  can  be  predicted  and  compensated  for  in  real  time,  provided  that  the 
beamsplitter  temperature  is  known  and  that  the  responsivities  ratio  is  well  characterized. 
Fortunately,  the  CATSI  system  already  contains  a  temperature  sensor  installed  on  the 
beamsplitter  mount.  The  responsivity  ratio  can  be  accurately  evaluated  by  application  of  the 
two-temperature  calibration  method.  Moreover,  if  the  temperature  dependence  of  the 
beamsplitter  spectral  responsivity  is  predefined  (from  proper  measurements)  then  the  full 
calibration  of  CATSI  spectra  can  be  performed  in  real  time  without  any  additional 
measurements  which  reduces  by  a  factor  of  3  the  number  of  operations  necessary  to  generate 
a  calibrated  spectrum. 

Finally,  from  the  previous  discussion  and  by  extension  of  eq.  32,  it  follows  that  the 
calibrated  spectrum  generated  in  real  time  with  the  CATSI  prototype,  for  the  single-beam 
operation  mode,  is  given  by 

L2  =  K~2(TBS)  +  Ll  '  rC2  (TBS  }  [1°9] 


or  for  the  differential  operation  mode 


5  L  = 


K2(Tbs) 


-  ^(Tbs) 


[110] 


where  it  is  convenient  to  recall  that  S  is  the  measured  raw  spectrum  (complex),  K2  (Tgs)  is 
the  complex  responsivity  and  re2  (Tgs)  is  the  instrument  residual,  both  calibration 
parameters  being  determined  by  the  beamsplitter  temperature.  Although  the  CATSI 
prototype  was  conceived  to  be  fully  consistent  with  the  proposed  Fourier  spectrometer 
configuration,  it  was  not  possible  to  reach  the  ultimate  performance  (re2  =  0)  due  to  the 

optical  asymmetry  between  the  two  beamsplitter  substrates  discussed  above.  This  technical 
imperfection  in  the  beamsplitter  manufacturing  can  be  easily  eliminated  by  a  better  control  of 
the  optical  characteristics  of  the  paired  substrates.  Unfortunately,  there  was  no  other  CATSI 
beamsplitter  (of  this  type)  available  to  verify  this  asymmetry  hypothesis.  However,  in  the  next 


UNCLASSIFIED 

42 


WAVENUMBER  (cm'  *) 


WAVENUMBER  (cm'  *) 

FIGURE  13  -  Spectral  responsivities  (a)  and  offsets  (b)  associated  with  each  of  the  two 
optical  channels  of  CATSI  system.  Fig.  13c  compares  the  single  port 
measurement  of  a  blackbody  at  15°C  (upper  curve)  and  the  resulting 
difference  when  the  two  ports  are  utilized. 
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Section,  we  show  that  the  modeling  of  the  instrument  residual  based  on  the  beamsplitter 
asymmetry  agrees  very  well  with  the  observed  data. 

5.3  Modeling  of  the  Instrument  Residual 

To  evaluate  the  instrument  residual  re2  from  the  substrates  dissymmetry  discussed  in 
Section  5.2  (eq.  108),  we  first  use  the  results  established  at  eqs.  65b  and  67b  where  it  was 
found  that  the  raw  spectrum  associated  to  the  self-emission  of  each  substrate  was  given  by 

Ki  (-1-,;vT?lTsl)  bs  [mi 

I  cl  I  si 


and 


(1  -tc2ts2) 

Tc2Ts2 


Bs 


[112] 


for  substrate- 1  (fixed  comer  reflector  side)  and  substrate-2  (moving  comer  reflector  side), 
respectively.  It  is  remembered  that  Tc],  Tc2  and  Tsi  ,  Ts2  represent  the  transmittances  of  the 
associated  AR  coatings  and  substrates.  Bs  is  the  plank  radiance  related  to  the  substrate 
temperature.  In  this  case,  if  input- 1  and  input-2  look  at  target  sources  of  radiance  Li  and  L2, 
respectively,  the  resulting  raw  spectrum  is  given  by 

S  =  E1  +  E2  [113] 


with 


Ei  =  Kj 


Li 


+ 


U  -TciTsi) 
TclTsi 


[114] 


e2  =  k2 


l2  + 


(1  -tc2ts2)b 

Tc2Ts2 


[115] 


To  obtain  these  two  equations  it  is  assumed  that  the  main  source  of  uncompensated  (or 
unsuppressed)  self-emission  comes  from  the  beamsplitter  substrates  alone,  i.e.  other  self¬ 
emissions  from  the  two  input  components  cancels  out.  Rewriting  eq.  113  to  obtain  the  raw 
spectrum  in  terms  of  the  instrument  residual  re2  (see  eq.  105  )  yields 
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S  =  K2  (5  L  +  re2) , 


[116] 


or 


S  =  K2 


8  L  +  (l 

. 


Kl  ~  TclTsl)  B  . 
K2  TciTs1 


(1  -Tc2Ts2)  e 
Tc2Ts2 


[117] 


Further  simplifications  occur  by  using  the  results  established  in  eqs.  59  and  60  where  it  was 
found  that  the  ratio  of  the  two  responsivities  is  given  by 


Kj  _  TciTsi 

K2  Tc2Ts2 


[118] 


which  after  insertion  into  eq.  1 17  yields 


re2  =  (!  -|^-)[Ll  -Bs]  [H9] 

or  alternately 

re2  =  (1  +  ii)[Ll'Bs]-  [120] 


As  a  validation  of  this  modeling,  Figs.  14  and  15  compare  the  measured  and  predicted 
instrument  residuals  (re2)  in  the  mid-IR  and  the  far-IR.  The  measured  residual  is  obtained 

from  the  standard  two-temperature  calibration  method.  The  predicted  residual  (model)  is 
calculated  using  eq.  120  (model).  Figure  14  compares  measurements  and  predictions  (real 
and  imaginary  parts)  obtained  when  the  input  source  Li  is  a  blackbody  at  a  temperature  of 
Tb=5.9°  C.  In  Figure  15,  the  input  source  Li  is  a  blackbody  at  a  temperature  of  Tb=49.5° 
C.  For  both  cases  the  beamsplitter  temperature  was  33°  C  (Bs).  As  predicted  by  eq.  120,  the 
behavior  of  the  residual  strongly  depends  on  the  temperature  difference  (through  Li-  Bs) 
between  the  source  and  the  beamsplitter.  For  instance,  in  Fig.  14a  a  negative  temperature 
difference  of  -27.1°  C  (i.e.  5.9°  -  33°)  yields  a  negative  real  part  and  a  positive  imaginary  part 
of  the  residual  in  the  region  from  550  to  650  cm'1.  In  Fig.  15a,  this  is  the  reverse  situation  in 
which  a  temperature  difference  of  +16.5°  C  ( i.e.  49.5°  -  33°  )  yields  a  positive  real  part  and  a 
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WAVENUMBER  (cm'  ’) 


b) 


WAVENUMBER  (cm  ') 


c) 


WAVENUMBER  (cm'  ’) 


FIGURE  14  -  Comparison  of  the  measured  and  predicted  (model:  eq.  120)  instrument 

residual  obtained  when  the  input  source  Lj  corresponds  to  a  blackbody  at 
Tr=5.9  C  (a,b)  correspond  to  the  far  IR  and  (c)  corresponds  to  the  mid  IR 
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FIGURE  15  - 


a) 


WAVENUMBER  (cm'  ) 


b) 


WAVENUMBER  (cm1) 


C) 


WAVENUMBER  (cm' ') 


Comparison  of  the  measured  and  predicted  (model:  eq.  120)  instrument 
residual  obtained  when  the  input  source  L\  corresponds  to  a  blackbody  at 
Tb=49.5  C  (a,b)  correspond  to  the  far  IR  and  (c)  corresponds  to  the  mid  IR 
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negative  imaginary  part  of  the  residual.  Overall  the  model  (eq.  120)  appears  quite  accurate  to 
predict  the  behavior  of  the  instrument  residual.  Many  other  comparisons  indicate  the  same 
level  of  accuracy.  The  slight  discrepancy  between  the  models  and  observed  in  the  800-1400 
cm~l  is  believed  to  be  due  to  the  experimental  errors  present  in  the  measurements. 

5.4  Impact  of  Polarization  on  Responsivity  Phases 

In  Section  4.3  it  has  been  found  that,  for  the  proposed  beamsplitter  configuration 
(thin  layer  of  air  squeezed  between  two  thick  substrates),  the  phase  difference  between  the  two 
responsivities,  Ki  and  K.2,  must  be  180°  (out  of  phase,  see  eqs.  59-60)  at  any  wavenumber 

and  for  any  dissymmetry  in  substrate  transmissions.  However,  the  experimental  results 
obtained  with  the  CATSI  system  do  not  entirely  agree  with  this.  As  shown  in  Fig.  16,  in  the 
spectral  region  of  strong  substrate  dissymmetries  between  530  and  680  cm~l,  phase 
differences  other  than  180°  appear.  Figure  16  shows  an  example  of  the  phases  for  Ki  and 
K2  and  the  corresponding  difference  as  a  function  of  the  wavenumber.  In  the  following,  this 
anomaly  is  interpreted  as  an  effect  due  to  the  two  components  of  polarization  (s  and  p) 
interacting  with  the  beamsplitter  substrates.  The  equations  that  give  the  raw  spectrum  for  a 
CATSI-type  double-beam  interferometer  are 

S  —  Ei  +  E2 

with 

Ei  =  Ki  (L]  +  SEini) 

E2=  K2(L2  +  SEin2)  [ 123 ] 

where  K],  K2  and  SEjni,  SEin2  are  defined  as  the  responsivities  and  the  self-emissions 
associated  with  input- 1  and  input-2,  respectively  and  Li  and  L2  are  the  source  radiance 
terms.  Note  that  the  self-emission  of  the  beamsplitter  air  layer  is  Ebs  =  0-  Because  we  want 

to  focus  primarily  on  the  effects  of  polarization  we  neglect  the  self-emission  terms  in  the 
discussion:  we  pose  that  SEjnl  =  0  and  SEin2  =  0.  In  this  case  and  according  to  eqs.  59-60 

we  have 


[121] 

[122] 


Ei  =  Kj  e*  VLj 


[124] 


E2  =  K2  ein  e'V  L2  . 


[125] 
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The  two  previous  equations,  which  do  not  take  the  polarization  into  account,  express  the 
fact  that  the  two  complex  responsivities  are  180°  out  of  phase.  The  equations  for  the  case  in 
which  the  two  polarizations  s  and  p  are  taken  into  account  are  obtained  by  extension  of  eqs. 
124  and  125,  yielding 


Ei  =  E*  +  E j  =  (  K®  e5 vs  +  K^tfVP) Li 
E2  =  E|  +  Ej  =  (  Ks2  e*  v,  +  KPe;vP)  e*L2 


[126] 

[127] 


Consequently,  the  effective  responsivity  modules  and  responsivity  phases,  seen  by  a  system 
that  measures  the  sum  of  the  two  components  (such  as  CATSI),  can  be  defined  by 


Ei  =  KfVvi^L! 

E2  =  K|ff  e1  (v2cff  +  ^)  L2 


[128] 

[129] 


where 

K®ff  =  [  (Ksj ) 2  +  (Kj) 2  +  2  KS!  (  cos(vs) cos(Vp)  +  sin(vs)  sin(vp) ) ]  2 


[130] 


WAVENUMBER  (cm-1) 


FIGURE  16  -  Typical  responsivity  phases  and  the  corresponding  difference  measured  for 
the  CATSI  system 
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tan  (\i^ff)  =  Ki  sin(^)  +  sin(^p) 

K*  cos(Vs)+  coslvp) 

and 


[131] 


K|ff  =  [  (K|) 2  +  (KP) 2  +  2  K|  (  cos(vs)cos(Vp)  +  sin(ys)  sin(yp) )  ]  2 


pn  faff)  _  Kt  sin(y.)  -  KP  sin(V 
K|  cos(vs)  +  RP  cos(Vp) 


[132] 


[133] 


From  the  previous  equations  (131  and  133)  ,  it  can  be  shown  that  the  two  effective 
responsivity  phases  vie^  and  V2e*^  can  differ  when  one  of  the  three  following  conditions 

happen: 

1  -  Kis  *  K2s  and  KjP  =  K2P  and  Vs  *  Vp  - 

2-  Kjs  =  K2s  and  K]P  *  K2P  and  Vs  *  Vp  > 

3-  Kjs  yt-  K2S  and  KiP  &  K2P  and  Vs  ^  Vp  • 

The  first  two  conditions  arise  when  there  is  a  dissymmetry  in  the  optical  transmission  of  the 

two  substrates.  The  third  condition  involves  the  beamsplitter  phase  jumps  (in  transmission) 
which  differ  in  polarization  when  there  is  absorption  in  the  substrates  and  their  associated  AR 
coatings:  this  is  a  characteristic  of  the  Fresnel  coefficients  for  absorbing  interfaces  (see 
Heavens,  Ref.  8).  This  yields  a  Responsivities  Phase  Difference  (RPD)  equal  to 


RPD  -  (v2ff  "  +  71 )  [134] 

rather  than  RPD  =  n  obtained  when  the  two  substrates  are  non-absorbing  or  optically 
identical.  In  the  beamsplitter  of  the  CATSI  system,  the  effects  of  polarization  are  observed  in 
the  530  and  680  cm'1  region  (Fig.  16)  due  the  combined  absorption  of  ZnSe  substrates  and 
thorium  fluoride  presents  inside  the  AR  coatings.  Above  700  cm"1,  the  beamsplitter  substrates 
and  the  associated  AR  coatings  exhibit  a  low  level  of  absorption  yielding  a  RPD  =  7t. 
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5.5  Interpretation  of  Channel  Spectrum  in  Instrument  Residual 


The  aim  of  this  section  is  to  show  that  the  apparent  channel  spectrum  observed  in  the 
CATSI  instrument  residual  re2  arises  from  a  mismatch  between  the  thickness  of  the  two  ZnSe 

substrates.  This  channel  spectrum  is  more  clearly  observed  in  the  spectral  region  above  1800 
cm"1  as  shown  in  Fig.  17.  To  establish  the  origin  of  this  phenomenon  it  is  convenient  to  use 
the  model  of  spectral  residual  discussed  in  Sections.  5.2  and  5.3.  In  this  context  the 
instrument  residual  is  given  by  (see  eq.  108) 


[135] 


Figure  18  exhibits  the  ratio  of  the  two  responsivities  module  I  Kq  /  K2  I  corresponding  to  the 

case  of  Fig.  17.  Equation  135  indicates  that  the  channel  spectrum  observed  in  the  instrument 
residual  re2  comes  from  a  periodic  mismatch  between  the  two  responsivities.  As  found  in  eqs. 

59  and  60  this  responsivities  ratio  is  directly  related  to  the  transmission  of  each  substrate  such 
that 


Kl  =  Isl 

k2  ts2 


[136] 


where  it  is  assumed  that  the  transmissions  of  the  two  AR  coatings  associated  to  each  substrate 
are  identical. 


To  understand  the  cause  of  this  periodic  mismatch  between  the  two  substrate 
transmissions  Tsi  and  TS2,  we  consider  the  interference  created  by  the  multiple  reflections 

inside  each  substrate.  Moreover,  to  simplify  the  problem  we  consider  only  the  two  main 
components  of  interference  (two  beam  interferometry),  as  illustrated  in  Fig.  19.  Taking  into 
account  this  interference  yields  for  the  substrate  transmissions 

r  i37i 

Tsl  =  Ti  [  1  +  Rar  Rs  cos  (27ivndi)  ]  1  J 

Ts2  =  T2  [  1  +  Rar  Rs  cos  (27tvnd2)  ]  [13§] 

where  d ] ,  d2  and  Tj,  T2  are  the  thickness  and  the  single  pass  transmission  of  each  substrate 
respectively,  and  n  is  the  refractive  index  of  the  substrate  material  (n  ~  2.5  for  ZnSe).  Rar  and 


RATIO  OF  RESPONSIVITIES 
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FIGURE  17  -  An  example  of  an  instrument  residual  spectrum  (real  part)  showing  the 
presence  of  channel  spectrum 


WAVENUMBER  (cm'1) 


FIGURE  18  -  Ratio  of  the  two  responsivities  I  Ki  /  K2  I  (module)  showing  the  presence  of 
channel  spectrum 
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Rs  are  the  inner  reflections  at  the  air-coating  interface  and  the  air-substrate  interface, 
respectively.  The  previous  equations  tell  that  each  substrate  generates  its  own  high  frequency 
channel  spectrum.  Ratioing  the  two  substrate  transmissions  gives 

Ki  _  Tsi  _  Ti  [  1  +  Rar  Rs  cos  (27tvndi)  ]  [139, 

K.2  Ts2  T2  [  1  +  Rar  Rs  COS  (27tVnd2)  ] 

which  can  be  approximated  by 

«  Il[  l  +RarRs(Cos  (2:tvnd2)  -  cos  (27tvndi) ) ]  [14°] 

K2  T2 

where  only  the  first  two  terms  of  a  Taylor  expansion  (denominator  of  eq.  139)  was  conserved. 
Simple  trigonometric  transformations  of  the  previous  equation  leads  to 

1  -  2  Rar  Rs  sin  [7tvn(d2  +  di)]  sin  [7tvn(d2  -  di)] )  ■  t141l 

K2  T2 


FIGURE  19  -  Two  beam  interferometry  applied  to  each  beamsplitter  substrate  of  CATSI 
and  the  related  parameters 
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The  first  sinus  term  creates  a  high  frequency  modulation  in  the  responsivities  ratio  spectrum. 
The  corresponding  interfringe  is  given  by 


A  v 


1 

n  (di+  d2) 


[142] 


which  leads  to  A  v  =  0.2  cm'l  for  the  ZnSe  substrates  used  in  CATSI  (di  =  d2  =  1  cm  and  n 
=  2.5).  This  high-frequency  modulation  can  not  be  seen  in  the  CATSI  spectrum  because  the 
resolution  of  the  instrument  is  limited  to  1  cm"l.  The  second  sinus  term  is  the  cause  of  the 
channel  spectrum  (beating).  By  checkings  Fig.  18  we  can  evaluate  the  interfringe  to  be 
approximately  200  cm‘l,  yielding  a  thickness  mismatch  equals  to 


(  d2  -  di  ) 


— 1 - =  -_U -  =20  pm  . 

n  (A  v)  2.5  (200  cm-l) 


[143] 


This  result  agrees  quite  well  with  the  25  pm  of  thickness  tolerance  given  by  the  manufacturer 
of  ZnSe  substrates.  The  natural  conclusion  of  this  section  is  that  there  are  two  ways  to 
eliminate  the  channel  spectrum:  First,  by  improving  the  thickness  matching  of  the  two 
substrate  (ideally  to  better  than  a  few  microns),  and  second  by  improving  the  AR  coating 
performance,  i.e.  making  Rar  as  small  as  possible,  which  minimizes  the  multiple-reflection 

contributions  as  seen  in  eq.  141. 

6.0  EXAMPLES  OF  SPECTRAL  MEASUREMENTS  WITH  CATSI 


To  verify  the  optical  suppression  efficiency  of  the  CATSI  instrument,  a  series  of 
laboratory  measurements  have  been  performed  in  the  differential-detection  mode.  For  the 
first  series,  the  two  telescopes  were  aimed  at  the  same  target  source,  a  flat  plate  blackbody 
(12"xl2")  approximately  3  meters  away  from  the  interferometer.  Figure  20a  shows 
examples  of  instrument  residuals  (500  to  1500  cm-*)  obtained  for  three  different  blackbody 
source  temperatures  10.6°  C,  33°  C  and  49.9°  C.  Also  shown  for  comparison  purposes  is  the 
measured  radiance  (single  beam  mode)  corresponding  to  a  source  temperature  of  10.6°  C  . 
As  expected,  above  750  cm-1  the  instrument  residual  is  quite  insensitive  to  the  source 
temperature  and  is  near  zero.  Actually,  for  all  three  measurements  the  instrument  residual  is 
generally  smaller  than  2  %  of  the  blackbody  radiance  at  10.6°  C.  Below  750  cm**,  the 
suppression  efficiency  is  degraded  due  to  the  dissymmetry  in  the  absorption  of  the  two  ZnSe 
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700  800  900  1000  1100  1200  1300 
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500  1000  1500  2000  2500  3000 

WAVENUMBER  (cm'  *) 


FIGURE  20  -  Three  examples  of  differential  measurements  with  CATSI:  for  three  different 
blackbody  source  temperatures  (a),  for  three  different  chemical  vapors  (b) 
and  for  a  smoke  stack  at  300  m  (c) 
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substrate  (and  AR  coatings)  composing  the  beamsplitter.  As  seen  in  Section  5.3,  this  strong 
instrument  residual  can  be  properly  predicted  by  eq.  120  (or  108). 

The  second  series  of  laboratory  measurements  has  been  performed  to  study  the 
evaporation  characteristics  of  several  chemical  simulants  (readily  available)  having  DR. 
absorption  bands  similar  to  some  of  the  potential  chemical  threats.  For  that,  a  small-size 
vapor  generator  has  been  built  to  test  experimentally  the  overall  detection  methodology.  The 
actual  vapor  generator  assembly  is  made  of  a  0.5-m  long,  7.5  cm  diameter  pipe  covering 
entirely  a  100  ml  becher  containing  the  desired  chemical  in  liquid  form.  The  evaporation  is 
efficiently  realized  by  ventilation  (no  heating)  using  a  small  fan  embedded  in  the  middle  of 
the  pipe.  This  assembly  produces  a  steady  flow  of  chemical  vapor  at  a  concentration 
determined  by  the  speed  of  the  fan.  Figure  20b  shows  superposed  on  the  same  graph  the 
results  of  three  distinct  measurements  performed  on  CCI4,  methanol  and  acetone  vapors 

using  the  dual-beam  interferometer  operated  in  differential  detection  mode.  For  that,  the  two 
telescopes  were  aiming  at  different  parts  of  the  same  flat  plate  blackbody  (12"xl2")  which 
acts  as  a  uniform  background  at  a  selected  temperature  of  44°  C.  The  experimental  set  up 
was  such  that  the  generated  vapor  was  probed  by  only  one  telescope  and  was  covered  entirely 
the  FOV  of  the  instrument.  The  vapor  temperature  was  equal  to  the  ambient  temperature 
(approximately  24°  C)  providing  a  negative  temperature  contrast  of  20°  between  the  target 
vapor  and  the  background.  Figure  20b  shows  that  in  these  conditions  the  three  spectral 
signatures  are  easily  and  clearly  identified.  Note  that  for  these  results  the  small  instrument 
residua]  (re2)  was  removed  for  the  sake  of  clarity.  The  spectra  were  recorded  at  a  resolution 

of  4  cm-  5  with  a  recording  integration  time  of  4  seconds  (30  scans). 

The  last  measured  example  (4  cm-1  resolution)  presented  in  Fig.  20c  represents  the 
radiance  differential  from  the  smoke  stack  of  the  DREV  heating  plant  located  approximately 
300  meters  from  the  receiver  site.  The  observation  is  made  at  an  elevation  angle  of  ~1  degree 
having  a  cold  clear  sky  as  the  background  and  an  ambient  temperature  of  -15°  C.  In  Fig. 
20b,  the  small  instrument  residual  (re2)  was  first  measured  independently  and  then  subtracted 

from  the  measured  signature.  In  this  case,  the  radiance  differential  exhibits  a  rather 
conventional  smoke  stack  signature  with  a  strong  blackbody  emission  coming  from  the  hot 
parts  of  the  stack.  In  addition,  the  so-called  red  and  blue  spikes  associated  with  hot  CO2 

emissions  (edges  of  the  2350  and  650  cm'*  bands)  are  clearly  observed  together  with  the 
strong  emission  lines  of  hot  water  vapor  in  the  1200  to  1900  cm'1  region.  Others  spectral 
features  in  the  800-1200  cm'1  would  require  a  deeper  analysis  to  be  associated  with  real 
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species  or  to  a  spectral  background  mismatch  between  the  two  adjacent  scenes  used  for  the 
differential  detection. 


7.0  SUMMARY  AND  CONCLUSIONS 

An  improved  configuration  of  FTIR  spectrometer  for  the  remote  detection  of  vapor 
IR  emissions  has  been  proposed.  The  originality  of  the  approach  lies  in  the  use  of  a  double¬ 
beam  interferometer  optimized  for  optical  subtraction.  This  double-beam  configuration  is 
particularly  suited  for  target  differential  detection  where  the  simultaneous  suppression  of 
background  and  self-instrument  contributions  is  done  optically  and  in  real  time.  Another 
advantage  is  that  the  proposed  FTIR  spectrometer  requires  only  the  monitoring  of  instrument 
temperatures  to  generate  a  radiometrically  calibrated  spectrum  in  real  time,  which  represents  a 
definite  advantage  over  usual  FTIR  instruments  . 

The  optimum  FTIR  spectrometer  configuration  has  been  quantitatively  defined  to 
fulfill  three  design  conditions,  i.e.  balanced  responsivities  (input- 1  vs  input-2),  balanced 
emissions  (input-1  vs  input-2),  and  beamsplitter  transparency.  To  achieve  these  requirements 
we  have  proposed  a  fully  symmetrical  double-beam  interferometer  with  a  special  beamsplitter 
assembly.  The  beamsplitter  is  made  of  an  optically  thin  layer  of  air  squeezed  between  two 
thick  ZnSe  substrates  properly  covered  with  anti  reflection  coatings  on  their  external  faces. 

To  verify  the  proposed  configuration,  a  prototype  referred  to  as  the  Compact 
Atmospheric  Interferometer  Sounder  (CATSI)  has  been  developed  according  to  the  above 
specifications,.  Results  obtained  with  CATSI  have  shown  that  the  instrument  residual  is  not 
exactly  zero  but  usually  smaller  than  2  %  in  the  region  from  700  to  3000  cm'l  (3.5-15  pm). 
However,  an  unexpected  strong  residual  has  been  found  in  the  region  from  500  to  700  cm' 1 
(15-20  pm).  After  analysis  it  appears  that  these  imperfections  comes  from  the  limited 
performance  of  the  manufactured  beamsplitter.  In  the  strong-residual  region  (15-  20  pm), 
these  imperfections  are  partly  due  to  non-identical  antireflection  coatings  present  on  the 
external  face  of  each  substrate  and  partly  due  to  a  difference  between  the  transmission 
associated  with  each  substrate.  Both  effects  introduce  a  dissymmetry  in  the  self-emission  and 
the  attenuation  associated  with  each  inputs.  A  simple  model  which  predicts  the  instrument 
residual  in  terms  of  the  substrate  characteristics  has  been  utilized  to  explain  the  observed 
residual.  In  the  2%-residual  region  (3.5-15  pm)  it  has  been  shown,  based  on  a  simple  model, 
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that  the  observed  channel  spectra  can  be  explained  by  a  thickness  mismatch  of  the  order  of 
20  pm  between  the  two  ZnSe  substrates  forming  the  beamsplitter. 

The  experimental  results  obtained  with  the  CATSI  prototype  have  shown  a 
background  suppression  factor  of  approximately  98%  in  the  spectral  region  from  3.5  to  13 
pm.  The  differential  detection  capabilities  of  CATSI  have  been  illustrated  in  two  ways.  First 
to  evaluate  a  smoke  stack  (heating  plan)  located  at  approximately  300  meters  from  the 
reception  site,  and  second  to  evaluate  the  evaporation  characteristics  of  three  different 
chemicals:  CCI4,  methanol  and  acetone  vapors.  Based  on  these  measurements,  the  sensitivity 

of  CATSI  to  acetone  concentrations  has  been  estimated.  Future  work  will  be  oriented  to  the 
development  of  the  differential  detection  method  for  the  remote  monitoring  of  target  vapors 
and  to  field  validations  with  the  CATSI  prototype. 
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